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Abstract 

We reduce the construction of a weak solution of the Cauchy prob- 
lem for the Navier-Stokes system on R 5 to the construction of a solu- 
tion to a stochastic problem. Namely, we construct diffusion processes 
which allow us to obtain a probabilistic representation of a weak (in 
distributional sense) solution to the Cauchy problem for the Navier- 
Stokes system on a small time interval. Strong solutions on a small 
time interval are constructed as well 
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Introduction 

The main purpose of this article is to construct both strong and weak 
solutions (in certain functional classes) of the Cauchy problem for 
the Navier-Stokes (N-S) system in i? 3 . To this end we consider a 
stochastic problem and show that the solution of the Cauchy problem 
for the Navier-Stokes system can be constructed via the solution of 
this stochastic problem. 

The approach we develop in this article is based on the theory 
of stochastic equations associated with nonlinear parabolic equations 
started by McKean [1] and Freidlin [2], [3] and generalized by Belopol- 
skaya and Dalecky [4] , [5] on one hand and on the theory of stochastic 
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flows due to Kunita [6] on the other hand. In our previous paper [7] 
we have constructed a stochastic process that allows us to prove the 
existence and uniqueness of a local in time classical (C 2 -smooth in the 
spatial variable) solution of the Cauchy problem for the Navier -Stokes 
system. In the present paper we construct a process which allows us to 
obtain construction of solutions of the both weak and strong Cauchy 
problem for this system. Later we plan to apply a similar approach 
for the Navier-Stokes equation for compressible fluids extending the 
results from [9], [10]. 

A close but different approach is the Euler-Lagrange approach to 
incompressible fluids which was developed by Constantin [11] and 
Constantin and Iyer [12]. Shortly, the main differences in these ap- 
proaches are the following: we use a probabilistic representation for 
the Euler pressure instead of the Leray projection and obtain differ- 
ent formulas for the stochastic representation of the velocity field. We 
discuss these differences with more details in the last section of the 
present work. 

Thestructure of the present article is as follows. In the first section 
we give some preliminary information concerning different analytical 
approaches to the Navier-Stokes system. Here we recall some common 
ways to eliminate the pressure and to obtain a closed equation for the 
velocity. 

The classical approaches here are based on the so called Leray 
(Leray-Hodge)-projection that is a projection of the space of square 
integrable vector fields to the space of divergence free square integrable 
vector fields. Applying such a projection to the velocity equation 
one can eliminate the pressure p and get the closed equation for the 
velocity u. This operator is used both in numerous analytical papers 
(see [15] for references) and in papers where the N-S system is studied 
from the probabilistic point of view [16], [17], [12]. Finally the pressure 
is reconstructed from the Poisson equation. 

One more possibility to eliminate the pressure appears when one 
considers the equation for the vorticity of the velocity field u and uses 
the Biot-Sawart law to obtain a closed system. From the probabilistic 
point of view this approach was investigated in [18]. 

In our previous paper [7] we do not use the Leray projection but in- 
stead we start with consideration of a system consisting of the original 
velocity equation and the Poisson equation for the pressure and con- 
struct their probabilistic counterpart. The probabilistic counterpart 
of the N-S system was presented in the form of a system of stochas- 
tic equations. Furthermore we prove the existence and uniqueness 
of a solution to this stochastic system and show that in this way we 
construct a unique classical (strong) solution of the Cauchy problem 
for the N-S system defined on a small time interval depending on the 
Cauchy data. 

In the present paper we also reduce the N-S system to the sys- 
tem of equations consisting of the original velocity equation and the 
Poisson equation for the pressure but then an associated stochastic 
problem considered here allows to construct a generalized (distribu- 
tional) solution to the Cauchy problem for the N-S system. The as- 
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sociated stochastic problem is studied in section 5. In sections 1-4 
we expose auxiliary results used in section 5. Namely, in section 1 we 
give analytical preliminaries and recall the notions of strong, weak and 
mild solutions to the Cauchy problem for the Navier-Stokes system. 
More detail can be found the recent book by Lemarie-Rieusset [15]. 
In section 2 we give a short review of probabilistic approaches to the 
investigation of the Navier-Stokes system [7], [16] -[18]. In section 3 
we study a probabilistic representation of the solution to the Poisson 
equation, while in section 4 we recall some principal fact of the Ku- 
nita theory of stochastic flows and apply the results from [19], [20] to 
construct a solution of the Cauchy problem for a nonlinear parabolic 
equation (see also [21]). Finally all these preliminary results are used 
to construct the probabilistic counterpart of the Navier-Stokes system, 
prove that there exists a unique local solution to the corresponding 
stochastic system and apply the results to construct both the strong 
and weak (and simultaneously mild) solutions to the Cauchy problem 
for the Navier-Stokes system. 



1 Preliminaries 

As it was mentioned in the introduction the main purpose of this 
article is to construct both strong and weak solutions (in certain func- 
tional classes) of the Cauchy problem for the Navier-Stokes system via 
diffusion processes. 

Consider the Cauchy problem for the Navier-Stokes system 

— + (u, V)u = uAu - Vp, u(0,x) = u (x), xeR 3 , (1.1) 

divu = 0. (1.2) 

Here u(t,x) G R 3 ,x £ R 3 ,t G [0, oo) is the velocity of the fluid at the 
position x at time t and v > is the viscosity coefficient and p(t, x) 
is a scalar field called the pressure which appears in the equation to 
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enforce the incompressibility condition (1.2). Later we set v = ^- for 
reasons to be explained below. 

By eliminating the pressure from (1.1), (1.2) one gets a nonlinear 
pseudo-differential equation which is to be solved. There exist different 
ways to do it and we consider now some of them. 

Given a vector field / let P/ be given by 

P/ = / - VA^V ■ /. (1.3) 

Here and below we denote by u ■ v the inner product of vectors u and 
v valued in R 3 . 

The map P called the Leray projection is a projection of the space 
L 2 (i? 3 ) = L 2 (R 3 ) 3 of square integrable vector fields to the space of 
divergence free vector fields and we discuss its properties below. A 
quite direct definition of P is connected with the Riesz transformation 
Rj. Recall that R& = 7 fc A which means that for / G L 2 we have 
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T{Rjf) = where F(f) = / is the Fourier transform of /. 

Then P is defined on L 2 (i? 3 ) as P = Id + R <g> R or 

3 

(?f)j = fj + J2 R j R kfk- 

k=i 

Since RkRj is a C alder on- Zygmund operator, P/ may be defined on 
many Banach spaces. 
Set 

3 

j(t,x) = V k UjV jUk = Tr[Vn] 2 (1.4) 
k,j=i 

and note that 7 can be presented as well in the form 
7 = V- V- it<g>u = ^ VfcVj(nfcUj). 

By computing the divergence of both parts of (1.1) and taking into 
account (1.2) we derive the equation 

-Ap(t,x)=j(t,x) (1.5) 

thus arriving at the Poisson equation. The formal solution of the 
Poisson equation is given by 

p = A~ x 7 = A _1 V • V • u(g>u (1.6) 

since divit = and finally we present Vp in the form 

Vp = VA _1 V- V -U&U. 

Substituting this expression for Vp into (1.1) we obtain the following 
Cauchy problem 

du 

— = vAu- PV • (it u), u(0) = u o . (1.7) 
ot 

There are a number of ways to define a notion of a solution for 
the Cauchy problem (1.7). We will appeal mainly to the Leray weak 
solution [13] or to the Kato mild solution [14]. 



1.1 Leray and Kato approaches to the solu- 
tions of the Navier-Stokes equations 

Let V = V{R?) = be the space of all infinitely differentiable 
functions on R 3 with compact support equipped with the Schwartz 
topology. Let T>' be the topological dual of T> and denote by (</>, ip) = 
J R 3 (p(x)ip(x)dx the natural coupling between <j) £ T> and ip £ V . If it 
will not lead to misunderstandings we will use the same notation for 
vector fields it and v as well, that is 

3 

(h,u) = / ^h k (x)u k (x)dx. 
Jr3 fc=i 
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We recall that a weak solution of the N-S system on [0, T] x R 3 
is a distribution vector field u(t,x) in (D'((0, T) x R 3 )) 3 where u is 
locally square integrable on (0,T) x R 3 , div u = and there exists 
p G X>'((0,T) x i? 3 ) such that 

— = vAu - V ■ (u <g> u) - Vp, limu(i) = u (1-8) 

holds. 

The Leray solution to the N-S equations is constructed through a 
limiting procedure from the solutions to the mollified N-S equations 

r §s = Mi* - v ■ ((« * g £ ) ® u) - vp, 

^V-u = 0, (1.9) 
I lim t ^o u(t) = u . 

Namely it is proved that there exists a function 

u £ G L oo ((0,oo),L 2 )nL 2 ((0,T),(H 1 )) 

such that (at least for a subsequence u £k ) strongly converging in 
(L 2 oc ((0,T) x R 3 )) 3 to u which satisfies (1.9). 

Here H 1 is the homogenous Sobolev space H 1 = {/ € S' : V/ G 
L 2 } with norm ||/|| Hl = ||V/|| L2 . 

On the other hand to construct the Kato solution means to con- 
struct a solution u to the following integral equation 

u(t) = e tA u - C e('- s ' A PV • (n ® u)(s)ds. (1.10) 
Jo 

Note that instead of looking for u(t, x) and p(t, x) one can prefer 
to look for their Fourier images u(t, A) = (2-7r)~2 J RS e~ lX ' x u(t,x)dx. 

The Leray and Kato approaches stated in terms of the Fourier 
transformations of the Navier-Stokes system can be described as fol- 
lows. 

Applying the Fourier transformation to the relation (1.7) written 
in the form 

(h,u(t)) = (h,u(0)) + f\h,Au(s)}- f\h,V -(it® «)(«)} 

JO JO 

we derive the relation 

(h, u) = (h, u ) - f\h, \X\ 2 u(s))ds- (1.11) 

JO 

— ^-r f I I Y X k h l (X'),ui(s,X)u k (s,X- \')d\d\'ds. 
(2tt)2 Jo JRzJR 3 k % 

Here u corresponds to the Fourier transformation of u. 

On the other hand if we are interested in the Kato mild solution 
of the N-S system then we may apply the Fourier transformation to 
(1.10) and derive the following equation 

X (t, A) = exp{-u\X\H} X (^ X)+ (1.12) 
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2-u\X\Ht-s) 



^(x(s)ox(s)) 



(X)ds 



for the function 

3 



x(t^) = l(lY |A| 2 n(t,A). 



Here 

Xi°X2(A) = -^3 / fi3 (xi(Ai)-e A )n(A)x2(A-A0 |A/| l||^ A/|2 , (1.13) 
eA = rx[ and 

n(A) X = X-e A (x-e A ), (1-14) 

Coming back to (1.7) we note that the Leray projection allows 
to eliminate the pressure p(t,x) from the Navier-Stokes system, to 
construct u and finally to look for p defined by the solution of the 
auxiliary Poisson equation. 

Another way to eliminate p(t,x) from the system (1.1), (1.2) is to 
consider the function v(t,x) = curl u(i, x) called the vorticity. Since 
curlVp(i, x) = one can derive a closed system for u and v. Namely 
for u and v we arrive at the system consisting of the equation 

dv 

— + (u ■ V)v = isAv + (v ■ V)u, (1.15) 
and the so called Biot-Savart law having the form 

Here the cross-product u x v is given by 

d e 2 e 3 
u x v = det \ u\ U2 us 

, vi v 2 v 3 . 

(u 2 V3 - u 3 v 2 )ei + (u s vi - uiw 3 )e 2 + {uiv 2 - u 2 vi)e s , 

where (ei, e 2 , e%) is the orthonormal basis in R 3 . 

Note that the term (v ■ V)u can be written as (Vu)v or even as 
V u v , where V u is the deformation tensor defined as the symmetric 
part of Vu 

V u = ^(Vu + Vu T ), 
since by direct computation we see that 

(V u)v - V u v = -(V u + V u T )v = 0. 

To be able to present the precise statements concerning the ex- 
istence and uniqueness of solutions to the N-S equations we have to 
introduce a number of functional spaces to be used in the sequel. 
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1.2 Functional spaces 

We describe here functional spaces which will be used in the sequel. 

Let V = V(R 3 ) be the space of all infinitely differentiable functions 
on R 3 with compact supports equipped with the Schwartz topology. 
Let V be the topological dual to V. The elements of V are called 
Schwartz distributions. 

The space of i? 3 -valued vector fields h with components hk G T> 
shall be denoted by D(i? 3 ) and D' shall denote the space dual to 
D(i? 3 ). 

Let L q (R 3 ) denote the Banach space of functions / which are ab- 
solutely integrable taken to the q-ih power with the norm ||/|| ? = 

(J R3 \f(x)\idx)l, 

Let Z denote the set of all integers, and suppose that k £ Z is 
positive and 1 < q < oo. Denote by W k ' q = W k,q (R 3 ) the set of all 
real functions h defined on R 3 such that h and all its distributional 
derivatives V a of order \a\ = J2 a j < k belong to L q {R?). It is a 
Banach space with norm 




(1.17) 



We denote the dual space of W k ' q by W k ' m where ^ + - = 1. El- 
ements of W~ k,q can be identified with Schwartz distributions. The 
space W~ k ' q is also a Banach space with norm 

U\\-k, q = sup \{<f>,h)\, 

\W\k, q <l 

where 

(cf), h) = I 4>(x)h(x)dx. 

The spaces W k,p for k G Z and p > 1 are called Sobolev spaces. If 
p = 2 we use the notation H k for the Hilbert spaces W k ' 2 . In a 
natural way one can define the spaces W k ' q , H fc of vector fields with 
components in W k ' p , and H k and so on. 
Set 

V = {v £ D : divv = 0} 

and let 

H = {closure of V in L 2 (ii 3 )}, V = {closure of V in H 1 }. 

(1.18) 

Let C%{R 3 ,R 3 ) denote the space of k-times differentiable fields 
with the norm 

Ibllc- = E Halloo 

\(3\<k 

and let C*' a (R 3 ,R 3 ) be the space of vector fields whose k-th deriva- 
tives are Holder continuous with exponent a, < a < 1 with the 
norm 

HslU*'" = WdWct + id]k+a 
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where 

_ \DPg(x)-DPg(y)\ 

[g\k+a = }^ SU P 1 _ i a • 

w=k x,yen? \ x V\ 

We denote by Lip(i? 3 ) the space of bounded Lipschitz continuous 
functions with the norm 

n n \g{x) - g{y)\ 

\\g\\Lip = sup — _ — . 

x,y£R? \ x V\ 

Spaces of integrable functions on the whole R 3 appear to be not 
satisfactory to construct a solution to the N-S equations and one has 
to consider spaces of locally integrable functions. 

Let / : R? — ► R 1 be a Lebesgue measurable function. A set of 
functions {/ : J K \f(x)\ p dx < 00} for all compact subsets K in R 3 is 
denoted by L p oc and called a space of locally integrable functions. Note 
that L}(R?) C L\ oc {R 3 ) . Although L P 0C (R 3 ) are not normed spaces 
they are readily topologized. Namely a sequence {u n } converges to u 
in L p oc (R 3 ) if {u n } — > u in L P (K) for each open K C G having compact 
closure in R 3 . Local spaces W^(R 3 ) can be defined to consist of 
functions belonging to W k ' p (K) for all compact K C R 3 . 

A local space Wj£(G) is defined as a space of functions belonging 
to W k ' p (G') for all G' C G with compact closure in G. A function / G 
w!°J?(G) with compact support will in fact belong to Wq' p (G). Also 
functions in W 1,P (G) which vanish continuously on the boundary dG 
will belong to Wq(G) since they can be approximated by functions 
with compact support. 

In the whole space R 3 and with p, q satisfying 1 < q < p < 00 de- 
note by M p a nonhomogenous Morrey space and by M p a homogenous 
Morrey space with norms given respectively by 



M p = if € L q loc : \\f\\ M v = sup supRp * \\f \\ l,( B (x , R)) < 00 
[ x eR 3 0<R 



(1.19) 



and 
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M q=\fe L loc : \\f\\M* = SU P SU P Rt 5 ll/IU«(B(x 0) fl)) < 00 
{ x £R3 0<R<l ) 

(1.20) 

where B(xo,R) is a closed ball of R 3 with center at xo and radius R. 

Respectively the integrable function is said to belong to M q (G) if 
there exists a constant C such that 

/ \f{x)\dx < Ci? 3(1 "i } (1.21) 
JGnB R 

for all balls Br. The norm in M q {G) is defined as the minimum of 
the constants C satisfying (1.21) 

A distribution u on (0, T) x R 3 is said to be uniformly locally square 
integrable if for all 99 G D((0, T) x R 3 ) 

sup / / I \\(f(t, x — xo)u(t, x)\\ 2 dxdt < 00. 

rnpRsio JR3 
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Equivalently u is uniformly locally square integrable if and only 
if for all t < *i G (0,T) the function £/"t , tl (x) = (jj 1 ||u(t, x)|| 2 dt)2 
belongs to the Morrey space L 2 uloc . In this case we write 

uloc 

L 2 ((t ,ti) xi? 3 )). 

For 1 < p < oo the Morrey space of uniformly locally integrable 
functions on R 3 is the Banach space L p uloc of Lebesgue measurable 
functions / on R 3 such that the norm ||/|| P)U /oc is finite, where 

||/IU* oc = sup (/ \f(x)\*>dx)*. 

xoGR 3 ■'Hx— xo||<1 

For t < ti, 1 < p,q < oo the space L v uloc x L\ ((t , ti) x R 3 )) is the 
Banach space of Lebesgue measurable functions / on (to, h) x R 3 such 
that the norm 

sup(/ (T |/(t,x)|«dt)?dx)p 

a= eff J||a;-a;o||<l "'to 

is finite. 

A C 00 function / on R 3 is called rapidly decreasing if 
lim |Zr/ , (x)|(l + ||x||) ri = 

holds for any multi-index a and any positive integer n. Let S = 
S(R 3 ) be the space of rapidly decreasing C°°— functions equipped 
with the Schwartz topology and S' be the topological dual of S. Since 
5 includes V, S' is a subset of V. The elements of S' are called 
tempered distributions. 

1.3 Weak, strong and mild solutions of the 
Navier-Stokes system 

Now we are ready to give more precise definitions and statements con- 
cerning the existence and uniqueness of solutions of the N-S equations. 

Definition 1.1. (Weak solutions)^ weak solution of the Navier- 
Stokes system on (0, T) x R 3 is a distribution vector field u(t,x), u G 
(D'((0,T) x R 3 ) d such that 

a) u is locally square integrable on (0, T) x R 3 , 

b) V-u = 0, 

c) there exists p £ T>'((0, T) x R 3 ) such that 

dtu = Au — V • (u (g) u) — Vp. 

The classical results concerning the existence of square integrable 
weak solutions are due to Leray [13]. 

Theorem 1.1. (Leray 's theorem) Let uq G (L 2 (R 3 )) 3 so that 
V • u = 0. Then there exists a weak solution u G L°°((0, oo), (L 2 ) d ) n 
L 2 ((0, oo), (H 1 ) 3 ) for the Navier -Stokes equation on (0, oo) x R 3 so 
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that limt-,o\\u(t) — rto 1 1 2 = 0. Moreover, the solution u satisfies the 
energy inequality 

\\u{t)\\ 2 + 2( f \\V®u\\ 2 dxds < lluolll, (1.22) 

JO J R? 

where 

d ,. d d 

lK*)ll! = E l\u k (t,x)\ 2 dx, v®« = EEl^l 2 - 
k=i jR3 k=ij=i 

Definition 1.2. (Mild solution) The Kato mild solution of 
(1-1), (1-2) is a solution of (1.7) constructed as a fixed point of the 
transform 

r-t 

v 1 ^ e tA u {x) - / e ( *~ s)A PV • (v ® v)(9, x)d6 = e tA u - B(v, v). 
Jo 

(1.23) 

Note that the right hand side of (1.7) 

e tA u {x)- e ( '~ s)A PV- (u®u)(d,x)d6 = $(t,x,u) (1.24) 
Jo 

is a nonlinear map in the corresponding space and the solution u is 
obtained by the iterative procedure 

u ° = e tA u , u n+l = e tA u - B(u n , u n ). (1.25) 

Hence to construct a mild solution to (1.1), (1-2) means to find 
a suitable functional space for which &(t,x,u) given by (1.20) is a 
contraction. 

To this end one has to find a subspace £t of l? ulocx h\{(^,T) x 
i? 3 ) so that the bilinear transformation B(u,v) of the form (1.11) is 
bounded as a map £t x St £t- Then one may consider the space 
E T C S' defined by / G E iff / G S' and (e tA f) 0< t<T G S T and prove 
the following result. 

Theorem 1.3. The Picard contraction principle. 

Let St C L^ oca .L| ([0, T) x i? 3 ) be such that the bilinear map B is 
bounded on St Then: 

(a) If u G St is a weak solution for the Navier-Stokes equation 
(1.1) (1-2) then the associated initial value belongs to Ey. 

(b) There exists a positive constant C such that for all uq G Et 
satisfying V • u = and ||e* A uo||£ T < 00 there exists a weak solution 
u G St of (1.1) (1.2) associated with the initial value uq 

ft 

u = e tA u - e ( *" s)A PV- (u®u)ds. (1.26) 
Jo 

The classical results assert that for sufficiently smooth initial data 
for example for uq in the Sobolev space H k , k > 5 + 1, there exists 
a short time strong unique solution to (1.1), (1-2). On the other 
hand Leray proved the existence of a global weak solution of finite 
energy, i.e. u G L 2 called the Leray-Hopf solution. Although the 
uniqueness and full regularity of this solution are still an open problem 
nevertheless one knows that if a strong solution exists then a weak 
solution coincides with it. 
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2 Probabilistic approaches to the so- 
lution of the N-S equations 

Along with above functional analytical approaches recently a number 
of probabilistic approaches to the problems of hydrodynamics was 
developed (see [18]- [17], [7]). In this section we give a short survey of 
several different probabilistic approaches. 

Let (f2, J 7 , P) be a complete probability space, w(t), B{t) be a cou- 
ple of independent Wiener processes valued in R? . 

Assume that (u(t,x),p(t,x)) is a unique strong solution to (1.1), 

2 

(1.2) or (to be more precise) to (1.1), (1.5) and set v = ^- to simplify 
notations in stochastic equations. Since in this case u is C 2 -smooth 
one can check that the stochastic equation 

df( T ) = -u(t-T,£(T))dT + adw(T), £(0) = x, (2.1) 

has a unique solution and the function 

v(t, x) = E[uoW)) ~ f Vp(* - r, £(r))dr] (2.2) 
Jo 

satisfies (1.1) and hence equals to u by the uniqueness of the strong 
solution to (1.1). The relation 

roc 

-2p(t,x) = Ej(t,x + B(6))d0 (2.3) 
Jo 

with 7 given by (1.4) allows us to verify that divn = 0. Thus (2. 2), (2. 3) 
give the probabilistic representation of the solution to the N-S system. 

If u(t,x) is a C 2 -smooth solution to (2.1)-(2.3), then Ito's formula 
yields that for v(9, x) = u(t — 6, x) 

v(t,ffl) = v(0,x) + - (v,V)v + ^-Av}(r,ar))dT+ 

f aVv(T,^T))dw(T). 
Je 

Then (2.2) and the relation u' T (t — r, x) = —v' t (t,x) yield 

C f du a 2 
Eu(0, at)) = v(t, x) - E[ ^ [— + (u, V)u - — An] (t - r, ^(r))+ 

Vp(* - r, Ur))]dr] + E[ f Vp(t - r, Ur))]dr}. 

Jo 

Finally it results from (2.2) that 

r* du a 2 
E[ J o + («, V)n - y An + Vp](t - r, 6(r))]dr] = 0. 

Since the latter equality holds for all t and x we deduce that (l.l)also 
holds. The relation u(0,x) = uq{x) immediately follows from (2.2). 

The system (2.1)-(2.3) is a closed system of equations and we can 
try to look for its solution. Then at a second step we will look for the 
connection between this solution and a solution of the N-S system. 
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This approach was realized in paper [7]. It appears that to prove 
the existence of smooth solutions to (2.1)-(2.3) we have to consider 
the stochastic representations for V-u and Vp along with this system. 

Using general results of diffusion process theory and in particular 
the Bismut-Elworthy formula [8] we note that heuristic differentiation 
of (2.1)- (2.3) leads to the relations 

V k Ui(t,x) = E[VjUoi(£(t))r]j k (t)- 
f — (V iP (t - r,£(r)) f m {9)d Wl (e))dT] (2.4) 

JO CT JO 

and 

dr]ik = -VjUi(t - r, £{T))r) jk (r)dT, rj ik (0) = S ik . (2.5) 

In addition by Bismut-Elworthy 's formula (integration by parts) 
we can derive from (1.17) the probabilistic representation for Vp(i, x) 

oo I 

-E[i(t,x + B(s))B(s)]ds. (2.6) 

The main results in [7] can be stated in the following way. 

Let V = («, V«), Vi = {V(t,x) : H^Hl^ < oo}, if § < r < 2 
and V 2 = {V(t,x) : \\V(t)\\ LLc < oo}, if r > 3, V = Vi n V 2 n C 1+Q 
and let M = C([0,T],V) denote the Banach space with the norm 

\\V\\ r ,a = SU Pte[0jT] [\\V(t)\\ V + [Vu(t)] a ]. 

Theorem 2.1.([7}) Assume that V(0) = V G V. Then there ex- 
ist a bounded interval [0, T] depending on Vq and a unique solution 
(£(t),u(t,x),p(t,x), 7](t),Vu(t, x)) to the system (2.1)-(2.5) belonging 
to M for each r G [0, T] . 

Theorem 2.2. ([7]) Assume that the conditions of theorem 2.1 hold 
and uo G C 2+a . Then there exists an interval [0, Ti], T\ < T, such 
that for all t G [0, Ti] there exists a unique solution to (1.1), (1-4) * n 
M c M where M = MPiC 2 and this solution is given by (2.2), (2.3). 

A close approach based on a similar diffusion process was devel- 
oped by Busnello, Flandoli, Romito [18], though their starting point 
was the system that governs the vorticity v = curln and velocity u. 
The corresponding probabilistic counterpart of (1.15), (1.16) can be 
presented in the form of the following stochastic system 

df (r) = - u (t - r, Z{r))dT + adw(t), f (a) = x, (2.7) 
and the following two relations 

v(t,x) = E[U(t,s)u (m), (2-8) 

roc I 

2u(t,x)= -E[v(t,x + B(6))xB(6)]dB, (2.9) 

Jo v 

where U(t,s) = exp(fgVu(t — r, ^(r)(ir) is a solution to the linear 
equation 

dUl' x = Vu(t - s, ^(s))f/*' x ds, Ul' x = Id. (2.10) 



2Vp(t,x) = - 

Jo 
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The main results in [18] are as follows. Denote by 



U a (T) = 

{uGC([0,T},Cl ) (R 3 ,R 3 ))nL^([0,T],C^ a (R 3 ,R 3 ))\dWu = 0} 
the Banach space endowed with the norm 



\\u\\u<* = sup ess < t <T\\u{t)\\ c i ia 

b 

and by 

V a ' q (T) = {v£ C([0,T],C b (R 3 ,R 3 )) nL°°([0,T},C^(R 3 ,R 3 ))} 

the Banach space endowed with the norm 

\\v ||v«.9 = sup ess < t <T\\v(t)\\L<ir\C£ ■ 

Theorem 2.3.([18]) Given p G [l,|),a G [0,1] and T > Zei 
Co G Cg{R 3 , R 3 ) n LP(i? 3 , i? 3 ) and set 

£o = lbo||c«nLp- 

Then there exists r G [0,T] depending only on Eq, such that there is a 
unique solution of (2.7)- (2. 10). The diffusion process £(t) pZays a ro/e 
o/ a Lagrangian path, vector field u belongs to U a , and vector field v 
belongs to V a ' p {r). In addition the deformation matrix U^ 1 satisfies 
(2.7). 

After these developments P.Constantin kindly attracted our at- 
tention to his papers [11] [12] where the Lagrangian approach was 
successfully applied to the investigation of the Navier-Stokes system. 
The presentation of this approach and the discussion of their similarity 
and difference will be given in the last section of the present article. 

A probabilistic representation of the solution to the Fourier trans- 
formed Navier-Stokes equation (1.11) was constructed by Le Jan and 
Sznitman [16]. 

To describe their approach recall the definition of the solution of 
the Fourier representation (FNS) of the Navier-Stokes system. 

First for a solution u(t, x) of (1.26) and its Fourier transform u(t, A) 
one can introduce a function x(t, A) defined on [0, T] x R 3 such that 

Xt(A) = ^(|)t|A| 2 u t (A), a.e. forte [0,71, 

and 

X*(A) • A = 0, xt(-A) = Xt(A). 
In addition, for Lebesgue a.e. A, Xt(A) solves the equation 

Xt(A) = exp(-v\\\ 2 t) Xo (\) + /* ^AlV^I 2 ^)!^ o X s(X)]ds, 

JO 2 

(2.11) 

where 



Xi o X2(A) = ~ J (xi(Ai) • e A )B(A)x 2 (A - A x )^- 



|A|dA 



! |A-A; 



2 ' 
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if the initial function xo is a measurable function and for a.e. A G 
i? 3 \{0} 



XO :R\{0}^C\ Xo (A)-A = 0, xo(-A) = xo(A) . 
Introducing the kernel X from i? 3 \0 to (i? 3 \0) 2 

| h(X 1 ,X 2 )K x (dX 1 ,dX 2 ) = ^ J fe (Ai,A-Ai) 2 

for /i > measurable on (i? 3 \0) 2 one gets 

Xi o X2 (A) = -i y (xi(Ai) • e x )U(X) X 2(.X2)K x (dX 1 ,dX 2 ). 

It turns out that K is a Markovian kernel with some remarkable 
features that allow to study existence and uniqueness problems for 
(2.11) with the help of a critical branching process on i? 3 \{0} called 
the stochastic cascade. Namely, LeJan and Sznitman have described 
a particle located in A such that after an exponentially holding time 
with parameter v\X\ 2 with equal probability the particle either dies 
or gives birth to two descendants, distributed according to K\. A 
representation formula for the solution of (2.11) is constructed as the 
expectation of the result of a certain operation performed along the 
branching tree generated by the stochastic cascade. 

One more probabilistic model was recently developed by M. Os- 
siander [17]. A binary branching process with jumps that corresponds 
to the formulation of solutions to N-S in physical space was con- 
structed in [17]. 

Once again the N-S system is reformulated incorporating incom- 
pressibility via the Leray projection P and then the Duhamel principle 
is applied to derive 



u = e utA uo 



- / f e -"(t-^)Ap V . (u®u){s)ds+ f e- u(t ~ s)A Vf{s)ds, 
Jo Jo 

(2.12) 

V • u = 0. (2.13) 

Let K(y, t) = (27rt) _ te _ ' 1 2^ be the transition density of the Brow- 
nian motion w(t) G i? 3 

(Py)ij = $ij ~ ( e y)i{ e y)j 
bi(y; u, v) = (u ■ e y )~P y v + (v ■ e y )P y u, 
62(2/; u, v) = h(y; u,v) + u ■ (I - 3e y e y )ve y . 
Then (2.12) can be rewritten in the form 

u(x,t) = u (x - y)K(y,2ut)dy+ 
Jrz 

[ [ {-j— ^-K(z, 2vs)bi(z: u(x — z, t — s), u(x — z, t — s)) + 
Jo Jrs 4 us 
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(^K(z, 2us) - -i-jj / K(z, 2us)dy) 

\Z\ 47r|z| 4 J{y.\y\<z} 

bi{z\ u(x — z,t— s), u(x — z,t — s)) + (K(z, 2vs)~P z — 

—^{I-Ze z e T z )( K{y 1 2us)dy)g(x-z,t-s)}dzds. (2.14) 

47r N J{v-\v\<z} 

Theorem 2.4. ([17]) Let h : R 3 -> [0, oo] and h : R 3 ^ [0, oo] 
iwtt /i locally integrable and h, h jointly satisfying 

/ h 2 (x — y)\y\~ 2 dy < h(x) and h(x — y)\y\~ l dy < h(x) 
for all x £ R 3 . If for all x € R 3 and t > 

h(x) 



(Airvt) 2 1 / uq(x — y)e dy\ < 



IXV- 
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and 

l!?(^)l<N 2 ^ 

then there exists a collection of probabilistic measures {P x :x£ R 3 } 
defined on a common measurable space (f2, J 7 ) and a measurable func- 
tion E : (0, oo) x — > R 3 such that 

P x {{uj : |E(t,w)| < /or a// 1 > 0) = 1 

/or which a weak solution u(t, x) to the N-S can be presented in the 
form 

u(x,t) = h(x) \ E(t,u)dP x (u) for allx e R 3 ,t> 0. 
Jn 

Furthermore the solution u is unique in the class 

{ve(S'(R 3 x(0,oc))) 3 :\v(t,x)\< 2 -^^} 

for all x£R 3 ,t>0. 

Our survey is still far from being exhaustive. As already men- 
tioned the discussion of the Euler-Lagrangian approach developed by 
Constantin and Iyer will be postponed to of the present paper the last 
section. 



3 A probabilistic representation of the 
solution to the Poisson equation 

Within the framework of the approach developed in this paper we 
intend to construct diffusion processes associated with the system (1.1) 
(1.5). First we will start with (1.5) and recall some results concerning 
the solution of the Poisson equation in an open domain G C R 3 . 
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First we recall that by the divergence theorem a C 2 (G) solution of 
— Ap = 7 satisfies the integral identity 



Vp • V(f> dx = — 70 dx 
G JG 



for all (f) ^ Cq{G). In the space Wq ,2 (G) which is the completion of 
Cq(G) under the inner product 



(p, 0) = / Vp • 



G 



the linear functional 



FU) = - ( j(f>dx 
JG 
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may be extended to a bounded linear functional on the space W ' (G). 
Hence by the Riesz theorem there exists an element p € Wq' 2 (G) 
satisfying (p,4>) = F{4>) f° r ah 4> £ Q}(G0- Then the existence of a 
generalized solution to the Dirichlet problem — Ap = 7 and p = 
on dG is readily established. The question of classical existence is 
accordingly transformed into the question of regularity of generalized 
solution under the appropriately smooth bounded conditions. 

We give in this section a brief summary of a probabilistic approach 
to the solution of the Poisson equation. We will try to give the prob- 
abilistic proofs of the necessary facts inasmuch as they are known. 
Proofs of similar statements can be found in [18]. The source for 
analytical results is the book by Gilbarg and Trudinger [22] . 

Consider the Poisson equation 

-Ap(x) = 7 (x) (3.1) 

where p and 7 are scalar integrable functions defined on G. A Newton 
potential with density 7 is defined by 

N - ,{x) = iIa¥^\C iv) '"'- <3 - 2) 

If 7 is regular and has a compact support then ./V7 is known to be a 
solution of the Poisson equation (2.1). 

To derive a probabilistic interpretation of the relation (3.2) we con- 
sider the generator A = ^A of a Wiener process B(t) € R 3 defined on 
a given probability space J 7 , P). It is well known that on the space 
Co(i? 3 ) of all continuous functions vanishing at infinity the Wiener 
process generates the strongly continuous semigroup 

T a (x) = £ 7 (x + B(t)), x e R 3 , t > 0, 7 g C {R 3 ). 

Given a function with a compact support in G we extend it to the 
whole space R 3 by zero. 

By a direct computation we can check that 

P OO P P OO ^ 

/ E[ 1 {x + B{t))]dt= l{x + y) - w e-±tM 2 dtdy = 

JO JB? JO (27Tt)2 
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/ iArT ] l{x + y)dy = 2N 1 . (3.3) 
Jr3 2ir\\y\\ 

To prove that p = 2Nj solves the Poisson equation 

— Ap = 7 

we need some additional regularity properties of Nj. 

Lemma 3.1. Lei 7 G L m (R 3 ) n L q (R 3 ) with l<m<3<q<oo. 
Then iV 7 G Cq^ 3 ) and 



II^tIIoo < <W|| 7 |U + 1| 7 | 



q)- 



Proof. First we note that for every l,m such that -j- + ^ = 1 by 
Holder inequality we have 

£| 7 (x + B(t))| = — i-y / |7(x + y)|e-^dy < (3.4) 
(2vrt)2 -/fl 3 



C4~^ + ^ll7llr < cu~£ 



sinC e -| + | = -|;. 



Finally we rewrite the left hand side of (3.1) as 

/"CO rl POO 

/ E[~/(x + B(t))]dt = E[ 7 {x + B(t))]dt+ E[j(x + B(t))]dt 
Jo Jo Ji 

and applying the estimate (3.4) for r = q and r = m we derive 

/ S[ 7 (ar + B(t))]rft<C(||7|U + 11711,), 
Jo 

with C = max(C m , C q ). 

By Sobolev embeddings it is known [22] that if 7 € L l {R 3 ) then 
iV7 G C(i? 3 ). To check that iV 7 G C (-R 3 ) we note that for any R > 
we can rewrite the left hand side of (3.3) as 

/•OO /'OO 

y E[ 7 (x + B(t))]dt = J o E[y(x + B(t))I llB{m>R ]dt+ (3.5) 

y o s[ 7 (x + B(t))i\ mm < R ]dt. 

Let us prove that the first term on the right hand side of (3.5) 
converges to uniformly in x as R — > 00 and the second term converges 
to as || 3; || — > 00 for each ii. For the first term we apply the estimate 
(3.4) to derive 

sup E[\-y(x + B(t))\I {mm>R} ] < 

c(\h\\ P + ll7ll ? )(*~^^{[i,oo)}W + r*/ {[0)1 ) } (i)) 

and 

sup E[\j(x+B(t))\I mm>R} ] < C*"f ||7|| ro ( I e-^dy)* 

xe-R 3 "'{||?/!I>-R} 
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as R — > oo. To obtain the estimate for the second term we apply (3.4) 
once again and obtain 

E[\ 7 (x + B(t))\I mtmR ]< 
Ct~^t{ h(y)\ p I { \\ y - x \\<R}dy)^I {[1 , oo)} (t)+ 

J R J 

cr%([ h(y)\ q i {ll y- xll < R} d y y«i m (t)), 

J R J 

that yields after the integration in time that the second term on the 
right hand side of (3.5) converges to 0, since 7 G L m (G) Pi L q (G) and 
is zero outside G. 

To study derivatives of A^7 we apply the Bismut-Elworthy-Li for- 
mula 

V Xi £7[ 7 (x + B(t))] = ±E[ 7 (x + B(t))Bi(t)] 

that holds for a regular 7. 

Lemma 3.2. Let 7 G L m (R 3 ) n L 9 ( 3 ) /or some l<m<§<3< 
g < 00. T/ien \7iV7 G C (i? 3 ) and /or eac/i x e R 3 

roc 1 

2V Xi N 7 (x)= -E[ 7 (x + B(t))Bi(t)]dt, i = 1,2,3. (3.6) 
Jo * 

Moreover 

IIva^tIU < c m9 (|| 7 |U + IItIU)- (3.7) 

Proof. By the Holder inequality 

1 C f 11 li 2 

-E\ 7 (x + B(t))B t (t)\dt = — h(x + y) yi \e-^dy< (3.8) 

°h\\ m t^ = c m h\\ m t- L i-%. 

Finally to give the sup estimate for the second derivative of iV~7 
one has to apply the Schauder estimates and the Bismut-Elworthy-Li 
formula. 

Let us recall two more useful results (see [22] theorem 4.5) con- 
cerning the Newton potential. 

Lemma 3.3.Let 7 G L q (R 3 ) n C£(R 3 ) with 1 < q < §. Then 
N-y€C} a (R 3 )nC {&), 

W N "f\\cf a m ^ c (hh*(R3) + ll7llc«(G)) 

and p = 2Nj is the unique solution of the Poisson equation 

— Ap = 7 

inC (R 3 )nC 2 (R 3 ). 

Theorem 3.4.Let G C^(i? 3 ),7 G C$(R 3 ) satisfy the Poisson 
equation ANj = 7 in R 3 . Then G R 3 and if B = Br(xo) is any 
ball containing the support of then 

||V 2 7V7|| 0iQ;S < CMo^-b, \\N 7 \\' hB < Ci? 2 ||7||o, B . (3.9) 
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In the sequel we will need as well L q type estimates for the Newtonian 
potential. 

Lemma 3.5. The operator N maps L q (R 3 ) into L q (R 3 ) and there 
exists a positive constant C such that 

\\Nj\\ Lqm < C||7l|L«(fl3) (3-10) 
Proof. By Holder inequality we have 

\N 7 \(x) = \ [ 7 (y)(r(x-y))i(r(x-y)) 1 -i^| < 
JG 

{[ \l{y)\ q T{x-y)dy}-«{[ T(x-y)dy} 1 ~^ < 
JG JG 

C{[ \~f(y)\ q r(x-y)dy}l 
JG 

Next we obtain by Fubini's theorem 

/ \Nj\P(x)dx< [ C p {[ \ 1 {y)\ p T{x-y)dy}dx = 

JR 3 JR 3 J R 3 

C p ! ! \ 1 {y)\ p T{x-y)dydx = C p f | 7 (y)| p ( / V{x-y)dx)dy< 

JG J R? JR 3 JR 3 



Ci 



[ h(y)\ p dy. 

JR. 3 



Note that all above results in this section are valid if we consider 
a bounded domain G C R 3 instead of R 3 . To get further regularity 
properties of the Newton potential we need more auxiliary results. 

Define the distribution ^ 7 (A) of the function 7 : G — ► R 1 by 

^ 7 (A) = |{xGG:| 7 (x)| >A}| (3.11) 

where |G| denotes the Lebesgue volume of the domain G. 

Lemma 3.6. Assume that 7 G L q (G) for some q > 0. Then 

f 7 (A) < A" 9 J |7(x)|dx, 
/ \Mx)\ q dx = p / A«-V 7 (A)dA. 

JG JO 
Proof. It is easy to check that 

/ \j{x)\ p dx > f \j{x)\ p dx > X p \{x : 7 (x) > A}| = X p uJ\). 
JG ->{7>A} 

If p = 1 we can apply the Fubini theorem to change the order of 
integration 

r r rh(x)\ roo r 

/ \j(x)\dx= / dtdx= / I{ xeGn (x)>x}dxd\ = 

•J G J G J J J G 

poo 

/ Vf{X)d\. 
Jo 
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For arbitrary q we have 



^ 7 «(A) = \{x : j(x) > \"}\ = ^ 7 (A<?) 

and hence 

p X q ~ l uJ\)d\ = v iq (\ q )d(\ q ) = / \j(x)\ q dx. 

JO JO JG 

Lemma 3.7. Let 7 G L 9 (G) /or some 1 < g < 00. T/ien iV7 G 
W^ q {G) and 

||V 2 iV 7 || L9(G) <C(g,G)||7||L 9 ( G ) (3.12) 
Moreover for q = 2 i/ie equality 



L 



V z Nj\\\x)dx = / 7 2 (x)cte (3.13) 
/ioWs. 

Proof. The proof of this fact is based on the Calderon-Zygmund 
technique of cube decomposition and estimates of the function ^ 7 (A) 
of the form (2.5). 

Let K be a cube in R 3 , 7 > integrable, and finally fix k > such 
that 

/ j(x)dx < K. 
\K\ Jk ~ 

Bisect K into 2 3 equal ( in volume) subcubes. Let Q be a set of those 
subcubes K for which f K "/(x)dx > k. For each of the remaining 
subcubes ( which do not belong to Q) we repeat the same procedure, 
that is bisect each one into 2 3 sub-cubes and add those smaller ones, 
where / is highly concentrated to Q. Now repeating the procedure 
again and again we obtain a partition of K. For any K in Q denote 
by K its immediate predecessor. Since K G Q, while K £ Q, we have 

1 f 1 f I K\ 1 f 

X< T1F\ l( x ) dx < /„ 7(a0<& = ^4-^7 /. 7 (x)da; < 2 3 A. 
Va- l-K-l Ja- |a| |k| 

Set F = \J KeQ K, J = K\F = r) K eQ.K c . Note that each point in J 
belongs to infinitely many nested cubes with bounded concentration of 
7 with diameters converging to zero, that is j^-y "y{x)dx < k, with 

\Ki\ — ► 0. By the Lebesgue theorem we deduce that p^-y fpr j(x)dx — ► 
7 a.e. with respect to the Lebesgue measure, that is 7 < ft a.e. on J. 
Then we have an average estimate on F and a point-wise estimate on 
J. 

At the second step we need the Marcinkiewicz interpolation theo- 
rem. 

Marcinkiewicz interpolation theorem.Zei 1 < q < r < 00 

and let T : L q {G) n L r (G) -► L«(G) n U\G) be a linear map. Suppose 
there exist constants C\,C2 such that V7 G L q (G) H L r (G) and for any 
A > 
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Then for any exponent m such that q < m < r the map T can be 
extended to a map from L m (G) to L m (G) and 

\\Ti\\l™-{G) < KCfC\ a \\^\\L m (G)- 

all 7 G L q (G)DL p (G) where ± = ^ + ±=^and the constant K depends 
only on m, q and r. 

At the end we define an operator T : L 2 (G) — > L 2 (G) by T7 = 
VjVj to obtain the necessary result. 

Theorem 3.8.(Calderon-Zygmund inequality) Let 7 G L P (G), 
1 < q < co . Then the Newton potential Nj = p G VF 2,9 (G) ; solves 
the Poisson equation Ap = 7 a. e. and 

l|V 2 p|U, (G ) < ChlU^G), (3-14) 
where C depends only on d and q. Furthermore , when q = 2 we have 

I \\V 2 N-f(x)\\ 2 dx = [ -f 2 (x)dx. 

JR 3 JG 

For the proof of the above interpolation theorem and theorem 3.8 
see, e.g., [22]. 

4 Probabilistic representations of 
weak solutions of parabolic equations 

In this section we adapt the results of the Kunita theory of stochas- 
tic flows acting on Schwartz distributions [19], [20] to the case under 
consideration. The considerations in this section are similar to [21]. 

Unlike the Kunita case we assume here that the coefficients of 
SDEs under consideration are at most C 1+a -smooth with < a < 1, 
but on the other hand it is enough for our present purpose to restrict 
ourself to nonsingular initial data for the Cauchy problem for parabolic 
equations and hence we consider stochastic flows in Sobolev spaces H k 
for k = 1, —1. 

It is more convenient for computational reasons to use sometimes 
the Stratonovich form of the Ito equation. Recall that a process £(t) 
having the Ito differential of the form 

d£(t) = [a(£(t)) + ±TrVa(at))<r(ti(t))}dt + o-m))dw 
has the Stratonovich differential of the form 

d S m = a(ffl)dt + a(£(t)) o dw. 

We say that condition C 4.1 holds if for all t G [0, T] the functions 
g(t) and a belongs respectively to Cl +a and C 2+a . 

Throughout this section we assume that C4.1 holds. We shall first 
give a brief review of the results which will be needed in the sequel. 

Consider a stochastic differential equation in the Stratonovich form 

d£(r) = ~g(t - r, e(r))dr - <t(£(t)) o dw(r), £(s) = x G R 3 , (4.1) 



21 



< s < t < t. Here g(t, x) G R 3 , a(t, x) G R 3 x i? 3 and G i? 3 is 
a Wiener process. 

Assuming that g(t) G C 1 (-R 3 ) and a(t) is a C 2 -smooth matrix we 
are in the framework of the Kunita theory [6] and know that there ex- 
ists a local C 1 -diffeomorphism of R 3 generated by the solution £ s ,x{t) 
of (4.1). 

Namely, by general results on the SDE theory the existence and 
uniqueness of the solution £ S) x(t) to (4.1) are granted for a C 1 - smooth 
bounded function g. Moreover, in this case, one can prove that the 
solution £f x (t) of (4.1) has a modification 0f (x,u;) such that for all 
ll> outside a null set N C £1 

1) </>§ t (x,lo) is continuous in (s,t,x), and differentiable in x; 

2) ^(^(x,^) = <% >T (x,u>), if < s < r < t; 

3) the mapping <jj% (ui) : i? 3 — > i? 3 is a C 1 - diffeomorphism in R 3 . 
The map </>§ T (u;) is called a stochastic flow of C 1 - diffeomorphisms 

in i? 3 . 

We will denote by (0§ )T ) _1 (w) = ifij. s (u}) the map inverse to the 
stochastic flow 0f (cj) and will write simply ip^ s (x) for 4>^ s (x, uo). We 
check a simple property of an inverse stochastic flow. 

Lemma 4.1. Consider the a-algebras 

= a {w(9) : 9 G [0, s]} ^ = <t{u>(t) - u)(n) : a < n < r < *} 

and a continuous bounded process m(s) adapted to Ff. Then the 
process f(s) = g(t — s) for s G [0, t] is Tf s adapted and for all a, (3 
such that 0<a<(3<twe have 

rP rt-P 
/ f{r)dw{T) = / g{s)dw{s). 



Proof. Note that since w(s) = w(t — s) — w(t) we have w(f3) — 
w(a) = w(t — (3) — w(t — a), that yields T'i_ s = Tf s . 
Now we consider a partition of the interval [0, t] 

{0 = t <h<...<t k < t k+1 <...<t N = t} 

such that \t k+ \ — tk\ —>■ as N — > oo. Set 6 k = t — t k for k = 1, . . . , N, 
then 

rP N 

/ f(s)dw(s)= lim ^2f{t k )[w(t k+1 )-w(t k )] = 

N 

iim a(t - e k )[w(t - e k+l ) - w(t - e k )\ = 

n ^°° k=l 

, t - a 

- lim V g(6k) [w(6 k +i ) - w(6 k )] = - g(s)dw(s). 

The main point of Kunita's theory is that the stochastic flow is 
a bijection and that the inverse stochastic flow satisfies a couple of 
SDEs which will be used for different purposes. One of these SDEs 
is given by the following lemma due to Malliavin ( see ([23], lemma 
5.2.2) or [5]). 
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Lemma 4.2. Let £, g (r,x,w) be a solution of the stochastic equation 
(4-1) with s = or equivalently of the SDE 

dH(r) = -g(t - r, £(t))o!t + m(£(r))dr - *(Z(t))<1w(t), £(*) = s 

(4.2) 

where m(x) = \TrV a(x)a(x) . Then, for every fixed T > we have 
£(t - 0, x, w) = |(0, £(t, x, «;),«;) 

/or ewery < < £, and x, a.s. (P w ). 

In what follows we need as well some generalizations of the ltd 
formula. The first one called the ltd -Wentzel formula reads as follows. 

Lemma 4.3. (ltd- Wentzel formula) Assume that the process 
£(t) G R 3 has a stochastic differential of the form 

<%(t) = g(t,£(t))dt + a($(t))dw(t) 

and the process f(t, x) £ R 3 has a stochastic differential 

df(t, x) = *(i, x)dt + $(t, x)dw{t) 

with the same Wiener process w(t). Let the vector field \I/(i, x) G 
R 3 and the operator field $(i, x) £ R 3 x R 3 be C 2 smooth in x and 
continuous in t. Then the process rj(t) = f(t,£(t)) has a stochastic 
differential 

df m (t, f (t)) = * m (t, f (t))dt + $mk(t, Z(t))dw k + Vi/ m (t, 

(4.3) 

2 ViVj/m(*, + Vi$ mfe (i, £{t))o-ik{Z(t))dt. 

Remark 4.4. Note that (4.3) can be rewritten in the Stratonovich 
form as follows 

4/m(*, £(*)) = *m(«, + *mk(*> + Vi/ m (t, £(*)) + 

(4.4) 

+V i * mfc (t,£(t))(7 i fc(£(t))di. 

We apply lemma 4.3 to check that the inverse flow ipf to the flow 
(f) g (0,t) (generated by the solution of the equation in (4.1)) can also 
be represented as a solution of the following stochastic equation 

#? ;0 (x) = V^ )t (^f i0 )-^(t, x)dt + V^^orM*) ° dw, (4.5) 

where (V0q t ) _1 is the inverse matrix of the Jacobian matrix V</>q t of 
the map <f>Q t . 

Namely, we have the following statement proved by Kunita (see [6] 
Theorem 4.2.2) in a slightly different context. 

Theorem 4.5. Let C4.1 hold and (f>Q t be the solution of the equa- 
tion (4-1)- Then the inverse flow [<^oJ _1 = Vfo satisfies (4-5). 

Proof. To verify the statement of the theorem note first that the 
Jacobian matrix n g (t) = V0g t solves the Cauchy problem for the 
stochastic equation 

dK(r) = -Vg{T^l T {y))K{T)dT-Va{4>l T {y))K{T)odw{T), k(0) = J. 

(4.6) 



23 



Then, consider the stochastic process 

G(y,t) = j\<fi iT (y)- 1 g(T,<ft tT (y))dT+ 

^ < V^ T (y)-V« T (y))o^(r), 
and compute ^o.tC^oO^))' where ipf has the stochastic differential 

#f )0 = dG(^f )0 (x),t). 
Set 4>oj(y) = 4> 9 (y,t). By the Ito-Wentzell formula we have 

MMoW) = x+^d S ^(^ (x),9) + [ V<PW ei0 (x),0)odtf 9i0 {x) = 
X ~ Jo 9{6, ttMoW' 9))d0 - Jo ff (^W,oW. )) dw (0)+ 

f v^(^ i0 (x), e)\y<fif>y$ t0 (x), e)]- l g (e, ^ 9 (^ (x)))de+ 

J 

jf* V^(^%(x), 0)[V<W* o (z), ^)]-V(^ i ,(^ i o(x), 0)) o d«,(0) = x. 

Hence, ^ t (ipf q(x)) = x and thus tpf is the inverse to <f>Q t . 

Remark 4.6. Recall that by lemma 4.2 the process ipt,o(x) = £(t) 
along with (4.5) satisfies the SDE 

^(0) = g(0, i(e)))dB + <7(|(e)) o cM)(0). (4.7) 

Denote by Jq 4 (o;) the Jacobian of the map 0Q t (u;). Given /i 6 
and / G one can define the composition of / with the stochastic 
flow ipfo(x) as a random variable with values in Hr x defined by the 
relation 

(s° t0 {u,),h) = (f, w^K>)>> heH ^ ( 4 - 8 ) 

for any t and W ^ N . Note that if / is a distribution of the form 
/ = f(x)dx where / is a continuous function then / o ip tj0 is just the 
composition of the function / with the map ipt,o(w) and 

/ f(i/>t,o(y,u))h(y)dy= f(x)h((/> 0t t(x,u))J 0j t(x,(j)dx 

by the formula of the change of variables. 

Remark 4.7. Consider the case of constant diffusion coefficient 
a(x) = a and assume that the drift possessed the property divg = 0. 
Then (4.8) has the form 

(S 9 tfl (u;),h) = {f,ho<j% tt (u>)), hen 1 , (4.9) 

since in this case Jq t [ui) = Id is the identity map. 
Consider a linear PDE 

f = L*/- 7 (t), /(0) = / , (4.10) 
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where 

L 9 = (g,V)+L , 
L f = -FijViVjf + mjVjf, 

and 

and g G C 1 is a given bounded smooth function and /o £ C 1 (or more 
generally /o G I?')- 

To construct a probabilistic representation of a weak solution to 
(4.8) in the case when the initial data /o is a C 1 function (or even a 
distribution /o G T>') we consider the stochastic process 

A(t) = /„- /*7(r)o^ dr (4.11) 

and define its composition with a stochastic flow ip 9 Q (x) solving (4.7). 
Recall that ip 9 (x) ^ s inverse to the stochastic flow 0o t (x) generated 
by the solution £(i) of (4.2). 

It is proved in [19] that the generalized solution of (4.8) is given 
by the generalized expectation of X(t) o ip 9 tQ {x). 

To define the generalized expectation we consider the Sobolev 
spaces W k,q or the weighted Sobolev spaces 5^ defined in section 
1 and check that EX(t) o ipf : o{ x ) 1S we ^ defined. 

Lemma 4.8. For each integer k and q > 1, T > there are exist 
positive constants Ck„ q ,T, c' k q T depending only on the flow ipt,o such 
that for any t G [0, T] 

E\\X(t) o ^o||* < c fe „ ? , r ||/ ||^ + 4„, )T f ||7(r)||^dr, (4.12) 

J 

for all f G W k ' q and 7 (t) G W k > q . 

If /o G W k ' q then by this lemma for any h G W~ k ' q there exists 

{S 0tt ,h) =E{X(t)oi> tfi ,h), 

and 5o,t can be considered as an element from W k ' q . This element 
will be called the generalized expectation of X(t) o^ tj0 and denoted by 
£?[A(t)o^ >0 ]. 

For A; = 1, q = 2 we consider (<Sf , /i) = o ipf , h) which is a 
continuous linear functional on H} and can be regarded as an element 
of Hr x . Set 

U t ,of = E[f o ^, ] (4.13) 

and call it the generalized expectation of / o i/j t . It is easy to see 
that Utfi is a linear map form Ti" 1 into itself. Moreover the family 
Ut, s f = E[f otp tjS ] possesses the evolution property Ut tT U T , s = Ut >s for 
any < s < t < T. It can be immediately deduced from the evolution 
properties of <j) 9 st and J 9 t . 

Finally we compute the infinitesimal operator of the evolution fam- 
ily Ut, s - To this end we need a version of the Ito formula. 
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Theorem 4.9. (The generalized ltd formula) Let f(t) £ Ti 1 be a 
continuous in t nonrandom function. Then, given stochastic flows 
00,4)^,0 generated by (4-1), (4-6) the following relations hold 

fit) o < t = /(0) + o $ >fl + L /W o $ e ]d0+ 

and 

fit) o = /(0) + j\^W)\ o < + L o[/(#) o ^, o p0+ (4.14) 

fv i [fie)o^]adw + f t V t [f(9)o^ eo ] g (8)d9. 
Jo 'Jo 

Here we understand the action of the operator Lq in the sense of gen- 
eralized functions. 

The proof of theorem 4.9 employs the classical ltd formula for 
C 2 - smooth functions f £ that approximate the C 1 function /, uses 
equations (4.2) and (4.6) for the flows <fo,t and iptfii respectively, and 
then justifies the passage to the limit under the integral sign in the 
integral identity. The details can be found in [19] for a much more 
general case. 

Let us come back to the parabolic equation (4.10) and set 7 = 0. 
We can show that the stochastic flow ip t ,o gives rise to an evolution 
family acting in spaces of of distributions and the function f(t) = 
E[fo ° ipt,o] is a weak solution of (4.10) with 7 = 0. 

Theorem 4.10. Assume that the coefficients of the stochastic 
equation (4-2) satisfy C 4.1 and ipf is generated by the solution of 
(4-6). Then, for any functions fo,g(t) G Ti 1 the relation 

f(t) = E[f otf fl ] 

defines the unique generalized solution to the problem (4-10) with 7 = 
0. The restriction ofUf to Ti 1 defines a strongly continuous family of 
evolution mappings acting on the space Ti 1 . The domain of definition 
of its infinitesimal operator A 9 (in a weak sense) contains the subspace 
Ti 1 and A 9 f = L 9 f for any u^Ti 1 . 

Proof. From the relation (4.13) and the properties of stochastic 
flows we deduce that the relation 

(US(t)f ,h) = (E[f o-fl fl ],h) 

defines a continuous linear functional on Ti 1 . Thus, we can treat 
U 9 (t)fo as an element from Ti 1 . It follows from the representation 
U^ fo(x) = E[foo ipf (x)] that U 9 (t) is a linear mapping from the 
space Ti 1 into itself. Note that the above definition of the family 
U% 8 through the integral identity allows to check that it possesses the 
evolution property 

(U t a T U? :S fo,h) = (U°J ,h). 
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Indeed, by the Markov property of the process V^o( x ) we deduce 
that 

(Ufjo, h) = (U 9 T U?J , h) = EKU&o, h o # jt J» t )] = 

E[(f , [vo^ s JJl T )\ v=ho<tja J=E(f , [ho^ T o^ T Pl T ]o^ s 4 :t ) = 

E[(fo, [hoct> 9 Stt ]Jl t )\ = {Uy,h). 
Now we apply the generalized ltd formula to obtain the relation 

E[fo ° fl, ] =fo + E[ f L 9 (f o ^ eo )d6}. 
Jo 

Note that in the latter expression each summand belongs to TC 1 . In 
addition, 

E[ /V(/o o h)d9\ = J\e[/ o ^ )0 ], (L»)*fc)d0 = 
j\L 9 (E[f o^ }),h)de, 

that yields 

£[/o°€,o] = /o+ /* Woo^ ])«w. 

JO 

In other words 

Uf,ofo = fo+ fuUlofodB. 
Jo 

As the result we get that /(£) = £[/o o ip t ] satisfies (4.8) and /(0) = 
/o- 

One can prove the corresponding result in the case j(t) 7^ in 
a similar way applying the above reasons to X(t) of the form (4.11) 
instead of /o- 

Theorem 4.11. Given tempered distributions /o and j(t) define 
\(t) by (4-11)- Then U(t) = E[X(t) o ipf ] defines the unique solution 
of equation (4-10) if ipffi satisfies (4-6) and </>o,t is its inverse. 

Proof. By the generalized Ito formula we get 

X(t) o Vf i0 = / - f l{r)dr + f Vi(A(r) o ^ )g(r)dr+ 
Jo Jo 

f Vi(A(r) o ^ )a(r)^(r) + /' Lg(A(r) o ^ )dr. 
jo jo 

As a consequence we get 

A(i) o Vf, = /o + T Vi(A(r) o # )<tcMt)+ 
Jo 

f L 9 {\{t)o^ )dr - [\(r)dT. (4.15) 

JO ' JO 
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Each term in (4.15) has a generalized expectation as an element of S'. 
The generalized expectation of the second term in the right hand side 
of (4.15) is equal to zero. For the third term we have 

(E \ f L 9 (X(r) o r T ,o)dr} ,h) = e\ f (X(t) o # >0 , [L 9 ]*h)dr 

J\e[X(t) o ^,oL [L 9 ]*h)dr = J\l 9 (E[X{t) o ^ ]), fc)dr. (4.16) 
Hence 

£[A(t) o Vf, ] = /o + T ^(^[A(r) o # ])dr - f 7 (r)dr. (4.17) 
Jo JO 

Differentiating each term with respect to t we check that U(t) = 
E[X(t) oipf fi ] satisfies (4.10). In addition lim t _> (U 9 (t), h) = (f ,h), 
that is lim t ^o(^ 9 (*) = fo and we proved that U 9 (t) solves the Cauchy 
problem (4.10). 

To prove the uniqueness of the solution to (4.8) suppose to the con- 
trary that there exist two solutions f(t) and f(t) to (4.10). Then the 
function u(t) = f(t) — f(t) satisfies = L 9 u{t) and \\m t ^Qu{t) = 0. 
Fix t and choose a function h(t, •) G V. Then there exists a solution 
h(r, x), < s < t < t, x e R 3 , to the Cauchy problem 

9h Q^ +[L 9 ]*h(r,x) =0, \im h(r, x) = h(t,x). 

If the coefficients a 9 and a 9 are C 1 -smooth , then there exists a unique 
classical solution to this Cauchy problem. As a result, 

«t), h) = j\^u(0), h{6))d0 + J\u(9), ^h(9)d6 = 

f \l 9 u{6), h{6))d9 - f\u(9), [L 9 }*h(9))d6 = 0. 
Jo jo 

5 A probabilistic approach to the 
Navier-Stokes system 

Let us come back to the Navier-Stokes system 



— + (u, V)u = — An - Vp, u(0,x) = u (x), x e R 3 (5.1) 
at 2 

-Ap = 7 , (5.2) 

with 7 defined by (1.4). 

Our main purpose in this section is to construct a diffusion pro- 
cess that allows us to obtain a probabilistic representation of a weak 
solution to (5.1), (5.2). To be more precise we intend to reduce the 
solution of this system to solution of a certain stochastic problem, to 
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solve it and then to verify that in this way we have constructed a weak 
solution of (5.1), (5.2). 

Let as above w(t),B(t) be standard R 3 - valued independent Wiener 
processes defined on a probability space (Q,J-,P). Given a bounded 
measurable function f(x) and a stochastic process £(t) we denote 
E StX f(£(t)) = Ef(£ StX (t)) a conditional expectation under the con- 
dition £(s) = x. 

In section 2 we recalled the probabilistic approach developed in our 
previous paper [7] that allows to construct a probabilistic representa- 
tion of a C 2 -smooth (classical) solution to (5.1)-(5.2) via the solution 
of the stochastic problem 

df( r ) =-u(t-T,£(T))dT + adw(T), (5.3) 

u(t,x)=E , x [u (at))+ f Vp(t - r,e(r))dr] (5.4) 

Jo 

roo poo 

-2p(t,x) = E[ i(t,x + B(t))dt]=E[ tr[Vu] 2 (t,x + B(t))dt]. 
Jo Jo 

(5.5) 

In this section we consider a similar stochastic system but now 
we choose to invert the time direction of the stochastic process itself 
rather then of the function u to obtain the possibility to reduce a 
construction of a generalized solution to the Navier-Stokes system to 
the construction of a solution of a stochastic problem. 

Our considerations will be based on the result of sections 3 and 
4. Note that since we consider the case where the diffusion coefficient 
a is constant the Ito form and the Stratonovich form of a stochastic 
equation coincide. 

Let as above w(t),B(t) be standard i? 3 -valued independent Wiener 
processes defined on a probability space (Q, J 7 , P). 

Let <f>o,t(,y) be a stochastic process satisfying the stochastic equa- 
tion 

d<h,t{v) = <f>o,t(y)) dt ~ <?dw(t), 0o,o = V 
and the stochastic process A(i) be of the form 

A(t) = «o- I* Vp{T,<h,r)dT. (5.6) 
Jo 

Consider the system 

dipt,e,i x ) = ~ u (9, ^t,e{x))d6 + adw{6), i/j t ,t{ x ) = x ' ( 5 - 7 ) 

u(t,x) = E[u (ip t ,o( x )) ~ f Vp(T,^tA x ))dr]. (5.8) 

Jo 

roc I 

-2Vp(t,x) = E[ -j{t,x + B(T))B{T)dT], (5.9) 

Jo T 

where 7 is given by (1.4) and prove the existence and uniqueness of 
its solution. 

To this end we apply the Picard principle to the solution of the 
stochastic system and construct a solution to (5.7)-(5.9) by the suc- 
cessive approximation technique. 
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Set 

u 1 (t,x)=u (x), ^ >0 (x) = x,p 1 (t,x) = (5.10) 

and consider a family of stochastic processes i^ k s (x) an d families of 
vector fields u k (t, x) and scalar functions p k (t, x) given by the following 
relations 

d^$ e = -u\e,$ fi )d0 + adw{O), $ >t = x, (5.11) 
u k+1 (t,x) = E[u (4 t0 (x)) - j\p k+l {T^l T (x))drl (5.12) 

/•oo 

-2p k+ \t,x)= E[ 7 k+ \t,x + B(T))]dT, (5.13) 
Jo 

where 

x ) = Tr[Vu k {t, x)Vu k+1 (t, x)}. (5.14) 

Note that for a fixed k the first stochastic equation (5.11) that de- 
termines the family of stochastic processes ip k (x) may be solved inde- 
pendently on equations (5.12)-(5.14). Then given the process ip k (x) 
and keeping in mind the properties of the function p k that satisfies 
the Poisson equation 

-Ap k+1 (t, x) = j k+1 (t, x), (5.15) 

one has to compute Vp k (t,x), u k+1 (t,x) by (5.12), (5.13). 

To investigate the convergence of the stochastic processes ip k (x) 
and functions u k (t,x),p k (t,x) defined above we need some auxiliary 
results concerning the behavior of solutions of stochastic equations. 

Let g € V be a given function. Consider the stochastic equation 

d^lg = -9(0) o -$ fi dB + adw(0), Vf» = x (5.16) 
and define vector fields u 9 (t,x) and Vp 9 (t,x) by 

U s(t,x) = E[u ^l (x)) - f Vp9(T^l T (x))dr}, (5.17) 

JO 

-2p9(t,x)= E[-f(t,x + B(r))]dT, (5.18) 
Jo 

~f 9 {t,x) =Tr[VgVu 9 ](t,x). (5.19) 
Recall that p 9 solves the Poisson equation 

-Ap 9 = Tr[VgVu 9 }. (5.20) 

To investigate the convergence of the stochastic processes ip k ,o( x ) and 
functions u k (t, x),p k (t, x) defined above we need some auxiliary results 
concerning the behavior of solutions of stochastic equations. Moreover 
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along with the system (5.16) - (5.18) we will need the system to de- 
scribe the process n k (r) = Vipt T ( x ) an d the functions Vu k (t,x) and 
Vp k (t,x). 

To derive the necessary apriori estimates we start with the consid- 
eration of a linearized system. 

Let g(t) be a given vector field. Consider the stochastic equation 

d^ = -g(e)o^ 9 d0 + adw(9), ^U x )= x ( 5 - 21 ) 
and define the vector fields u 9 (t,x) and Vp 9 {t,x) by 

u 9 (t,x) = E[u (r t ,o(x)) - [\p9(T^l T (x))dT], (5.22) 

J 

roo 

-2p 9 (t,x) = E[ 7 9 (t,x + B{T))]dT, (5.23) 
Jo 

~f 9 {t,x) =Tr[VgVu 9 ](t,x). (5.24) 
Finally we derive from (5.23) the relation 

roo I 

-2Vp 9 (t,x)= E[-^ 9 (t,x + B(T))B{T)]dT (5.25) 

JO T 

by applying the Bismut - Elworthy - Li formula first checking the 
conditions that validate such an application are satisfied. Below we 
will need some estimates of a solution of the Poisson equation from 
section 3. For convenience of references we formulate them in the 
following statement. 
Lemma 5.1. 

1. Let j 9 G L q (R 3 ) n L m (R 3 ) for some l<g<|<3<m<oo. 
Then 

||Vp 9 ||oo<C, m (|| 7 9 ||, + ||7 S IU) 
IIViV^lloo^Cdl^ll^ + Ma)- 

2. Let 7 9 G L r (R 3 ) for 1 < r < oo. Then p 9 G W^ r (R 3 ) and the 
Calderon- Zygmund inequality 

\\ViVjP g \\r,loc<Cl\\l 9 \\r,loc 

holds. 

Let Lip be the subspace of the space C(R 1 x R 3 ,R 3 ) of continuous 
( in t G [0, T], x G R 3 ), bounded functions which consists of Lipschitz- 
continuous (in x) functions g such that 

\\g(t,x)-g{t,y)\\<L g (t)\\x-y\\, t € [0,7] x,y£R 3 , 

where || • || is the norm in R 3 . 
Condition C 5.1 

Let g(t, x) G -R 3 be a vector field defined on [0, T] x R 3 that belongs 
to C 1,a (-R 3 , R 3 ),0 < a < 1 for a fixed t G [0, T] and satisfies the 
following estimates: 
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1- ^ N g (t) f° r some Q to be specified below, ||g(t)||oo < 

Zoc 

and 

||s(t, y)|| < L g (t)\\x-y\\, \\Vg(t, x)-Vg{t, y)\\ < L g {t)\\x-y\\. 

2. ||Vflr(*)||oo < Kl(t), ||V^(t)|| r , i0C < Nfo), where K g (t),L g (t), 
N g (t) and K g (t), L g (t), N g (t) are positive functions bounded on an in- 
terval [0, T], r = m and r = gforl<g<|<3<m<cx). 

Set ip(r) = ^t, T (x) and consider the stochastic equation 

ip( T ) = x-J g(Ti,if)(Ti))dTi +J odw{Ti), (5.26) 

with 0<r<t<Tfora certain constant T. If we are interested in 
the particular dependence of the process ip(r) on the parameters t, x 
and g, we write ip(r) = ipf,x( T )- 

Lemma 5.2. Assume that C 5.1 holds. Then there exists a unique 
solution ^|(r) of (5.21) that satisfies the following estimates: 

E\\r x (r)\\ 2 < 3[|M| 2 + a\t -T) + (t-r) j'iK^dn], (5.27) 

E\\r x {r)-r y {r)\\ < \\x-y\\ef>^ M , (5.28) 
E\mr)-^(r)\\ < f*\\g(T 1 )-g 1 (T 1 )\\ oo dT 1 ef' L °W d0 . (5.29) 

Proof. The proof of the estimates of this lemma is standard. We 
only show the proof of (5.28). In view of C 5.1 we have 

E\m(r) ~ r y (r)\\ < \\x - y\\ + J* L ff (ri)||^(n) - ^(n)!!^, 

where < r < t < T with some constant T to be chosen later. Finally, 
by Gronwall's lemma, we get 

EW x {r)-r y {r)\\<\\x-y\\Jt L ^ de . 

Along with the equations (5.21)-(5.23) we will need below the equa- 
tions for the mean square derivative n(t) = Vipt.o(x) of the diffusion 
process i/jt,o(x) that satisfies (5.21) and the gradient v(t, x) = Vu(t, x) 
of the function u(t,x) of the form (5.22). 

Lemma 5.3 Assume that C 5.1 holds. Then the process r/ 9 (r) = 
VV^r satisfies the stochastic equation 

d v °(r) = -Vg(r, ^f >T )^(r)dr, rf(t) = I, (5.30) 

where I is the identity map. Furthermore the process t) 9 {t) possesses 
the following properties. 

The determinant detr/(T) is equal to 1, i. e. 

detlf{T) = J t , r = l, 

and the estimate 

\\V 9 (r)\\<eIr K » e (5.31) 
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holds. 

In addition the following integration by parts formula is valid 

[ f(r x (r))dx=[ f(x)dx, feL\R 3 ). (5.32) 

Proof. Under the above assumptions the first statement immedi- 
ately follows from the results of the stochastic differential equation 
theory. By a direct computation one can check that Jj )T satisfies the 
linear equation 

dJ tyT = -div g(ipl T )Jt, T dT, J t ,t = I 

and since div g = we get the second statement that yields the in- 
tegration by parts formula (5.32). Finally (5.31) is deduced from the 
inequality 

E\\ V (t)\\<1 + j^Kl{0)EUe)\\dB 

by the Gronwall lemma. 

In the sequel we denote by rj x,9 (t) the solution of the equation 

dif> B (T) = Vg(t, ^ T (x))v x ' 9 (r)dr, r?^(0) = I 

if we will be interested in the properties of the process rf' 9 (t). One 
can easily check that 

W' g {T)-if>«{T)\\ < J* \\Vg(e,^ e (x)) - Vg(e,^ e (y))\\d9efr K » e 

< J\ l g (9)U t Ax))-Ae(y))\\de 

and by (5.28) we have 

\\v X ' 9 (r)-v v ' g (r)\\<C(r)\\x-y\\ 

where C(r) is a bounded function over a certain interval [0, Ti] de- 
pending on g. 

Let us state conditions on initial data no of the N-S system. 
We say that C 5.2 holds if for < a < 1 the initial vector field 
uq G C 1,a satisfies the following estimates 

1 1^0 1 |oo 

<K , ||V«o||oo<-Ko. \\uo\\r,loc < M , \\Vu \\ rMc < Mq 1 

with r to be specified below and let Lq,Lq be Lipschitz constants for 
the functions uq and Vuq respectively. 

Lemma 5.4. Assume that g(t,x) satisfies C 5.1 and uq satisfies 
C 5.2 with r = q and r = m for l<g<|<3<m<oo. Then the 
vector field u 9 (t,x) given by 

vP{t,x) = £MV t yO)) - f\p°(T^l T (x))dT] (5.33) 
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satisfies the following estimate 

IK(t)||oc < K + rC 9ro [||V(/(r)V« fl (r)||, ) , oc +||V5(r)V« fl (r)|| m) , oc ]dr. 
Jo 

(5-34) 

The proof of the estimate can easily be obtained by direct compu- 
tation from (5.33) using the estimates of the Newton potential given 
in lemma 5.1. 

Lemma 5.5. Assume that conditions of lemma 5.4 hold. Then 
given the function u 9 (t,x) of the form (5.33) the function Vu 9 (t,x) 
admits a representation of the form 

Vu°(t,x) = E\VuoW% fi (x))rf*{t)- 
f -77 —,Vp°(T, Yf iT (*)) f\^(9)dw(8)dr} (5.35) 

JO a(t — T) J T 

and the estimate 

\\VuV{t)\\ 00 <Jl K l^Kl+ 
f C qm — ^e/>«^4(r)[||V^(r)||^ oc + \\Vu^r)\\ m , loc ]dT 

JO O-yJt - T y 

(5.36) 

holds for l<(/<!<3<m<oo. 

Proof. To derive (5.35) we compute directly the gradient of the 
first term in (5.34) and apply the Bismut-Elworthy -Li formula [8] to 
compute the gradient of the second term in this relation. To verify 
the estimate (5.36) we use the above estimates for the process n(t) 
and the estimates of the Newton potential derivative from lemma 5.1. 
Then we obtain 

||V^(t)||oo<e/o^(^Ki + 
fc m , q —. ^efr K ^ )dd [\\Vg(T)Vu°(T)\\ m , loc + (5.37) 

JO 0~\Jt — T 

||V^(r)Vu»(r)||, )Joc ]dr] < eSo K l^ M K^+ 
fc^^^ef^^KlirMVu^r)^ + || W>(t)|L,^t]. 

JO 0~\/t — T 

Now we have to derive the estimate for the function || Vu(t)\\ r j oc . 

Lemma 5.6. Assume that the conditions of lemma 2.4 hold. 
Then for 1 < r < oo the function u 9 (t,x) given by (5.33) satisfies the 
estimate 

||V« ff (t)|| r , Joc < efo K * {e)de \\Vu \\r,ioc+ (5.38) 
C qm f do K M M Kl( T )\\ Vu°(r)\\ r , loc , dr 

with a constant C depending on r and a certain constant T to be 
specified later. 
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Proof. Recall that along with (5.35) Vu 9 (t,x) admits the repre- 
sentation 

Vu 9 (t,x) = E[Vu {ipf (x))ri x ' 9 (t) " C y 2 ^{r^ 9 tT {x))r] x ^{T)dT]. 

Jo 

To derive the estimate for ||Vtt(i)||£ Zoc = J G \\Vu(t, x)\\ r dx (where G 
is an arbitrary compact in R 3 ) by the triangle inequality we get 

\\^U 9 (t)\\r,loc < ai + « 2 , 

where 

ai = (| G ^[||Vn (^o(x))^(i)ir]^) F , 

a 2 = (j^ jf || W(r,< T (x)))7f ^(rJirdrdx) " . 

To estimate a\ we apply the Holder inequality and recall that 
ipt,r(x) preserves the volume. As a result we have 

ai < (jr,(S[||Vtio(<o(^))ll 2 ]^[ll^(*)ll 2 ]) 5 ^))' < 

llVtioll^^ 1 ^. 

To estimate «2 we apply the Calderon-Zygmund inequality and 
the above property of ipt jT (x) to obtain 

a\<C r f* eIo K l^ de Kl(r) [ || Vn(r, xWdxdr. 
Jo y Jg 

Combining the above estimates for a± and 0:2 we obtain the required 
estimate 

||V^(t)|| ri , oc <e/o^ 1 W« w [||V«o||r,Joc+ 
C-rj^Vo'^W^WIIV^WII^dr]. 

Theorem 5.7. Assume that conditions C 5.1 one? C 5.2 /ioW. 
T/ien i/iere exists an interval Ai = [0,Ti] and functions a(t), (3(t), 
k bounded for t € Ai, suc/i that, if for all t € Ai, ||g(t)||oo < K (^) 
and ||V#(t)||oo < a(t), \\Vg(t)\\ r < (3 r (t) then the function \\Vu 9 (t,x)\\ 
(where u 9 (t,x) is given by (5.21)) satisfies the estimates 

IKWHoo < «(t), ||V^(t)||L < a(t), ||Vu»(t)|£ Joc < r (t) 

(5.39) 

/or r = q and r = m and l<m<|<3<g<cx). 

Proof. Analyzing the above estimates for the functions u 9 (t,x) 
and Vu 9 (t,x) we get the following estimates 

HW^Hoo < efi K i ( - 0)d0 K%+ (5.40) 

^C^e^^W^i^i^fUV^C^)!!^ + ||V^(^)|| Woc ]dx], 
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||Vu(t)|| r)Joc < Jo K ^ )de [\\Vu \\ r , loc + (5.41) 
C r f elo K l^ de Kl(T)\\Vu3(r)U oc dT}. 
To derive the required estimates consider the integral equations 

a{s) = e /> ( ^K 1 + C , gm ^e/> (e)de a(r)K(r) + n m (r)]dT, (5.42) 

n r {s) = e^ Q{e)de ||Vu || r + C r j\$s a{e)de n r {T)a{T)dT 

for r = q and r = m and C qm = max(C q ,C m ). Finally we consider 
the equation 

(3( S ) = Jl a( ^ de (3 + C\ m J* J T s a ^ de a (r)(3(r)dT, (5.43) 

where f3(r) = n q (r) + n m (r), and 

||Vu ||^ oc + ||Vu ||m,; oc = n q (0) + n m (0) = /Jo- 
Next instead of the above system of integral equations we consider the 
system of ODEs 

^ = - a \s) - C qm a(s)P(s), a(t) = Kq, (5.44) 

^ = _ a ( s)/3 ( s )_C 9 1 m a( S )/3( S ), 0(t) = P<>- M5) 

By classical results of the ODE theory we know that there exists an 
interval [0, T\] depending on Kq, Nq and C, C qm such that the system 
(5.44), (5.45) has a bounded solution defined on this interval. 

To prove the convergence of functions u k (t,x),Vu k (t,x) we need 
one more auxiliary estimate. Actually, we have proved that u k (t) G 
Lip with the Lipschitz constant independent of k. It remains to prove 
that X7u k (t) have the same property. 

Lemma 5.8. Assume that C 5.1 and C 5.2 hold. Then the 
function Vu 9 (t) satisfies the estimate 

\\Vu°{t,x)-Vu°{t,y)\\<N?(t)\\x-y\\ a i/t e [0,31] 

for any x,y G G where G is a compact in B? and the positive function 
Nf(t) depending on parameters in conditions C 5.1 and C 5.2 is 
bounded over the interval [0,Ti] defined in theorem 5.7. 

Proof. Applying the integration by parts Bismut - Elworthy - Li 
formula to (5.33) we deduce the following expression for the gradient 
of the function u(t, x) 

Vu 9 (t,x) = E[Vuo(ipl (x))rj x ' 9 (t)- (5.46) 
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It results from (5.46) that 

\\Vu 9 (t,x) -Vu 9 (t,y)\\ < Ki + K 2 + K 3 + K 4 , 

where 

«:! = ^[||V^o« (^)) - V^ «o(y))llll^(*)ll], 
K2 = E[\\Vu (^ (y))\\\\ri x > 9 (0)-r,y>9( 

K3= f -tr^ E [W{r^ g tA x ))- 



/o cr(t - r) 
Vpf(r,^ r (y))|||| j\*> 9 {9)dw{9)\\]di 



r* 1 
k 4 = / — — i£ 

o 



' r 

ft 



l|V^(r,^ r (y))|| I \^ 9 {9)-^ 9 {9)]dw{9)\\ 



dr. 



a(t - t 

One can easily check using the estimates stated in lemmas 5.3 - 
5.5 that under conditions C 5.1, C 5.2 

Kl < LlE\\r t , (x) -Mv)\\efi Klme < \\x-y\\LlJ>^ 

and 

^2 < KlE\\^(t) - rf> 9 (t)\\ < \\x - y\\ K l g (9)J'l K ^ 0l)dei d9 . 

To derive the estimates for K3 and K4 we recall (see lemma 5.1) 
that the solution of the Poisson equation —Ap 9 = j 9 satisfies the 
estimates ||ViVjP»||oo < C(||7 fl ll« + [j 9 ] a , G ), ||V,V^|| r < || 7 »|| r an d 
1 1 Vp 9 1 |oo < C gm (||7 9 || g + 1 1 7 s 1 1 m ). Hence we obtain the inequalities 

*3 < f —?=(EUl T (x) - ^ r (y)|| 2 )^||7 s (r)|| 9 + [l 9 {r)] a , G ) 
Jo cryt — t 



and 



ay/t - t 

i:^dr< f t ^ ={ \\ x -y\\Llei: L ^( m r)\\ q+ 

JO ay/t — T 

b 9 {r)] a , G )Jl K ^dr 

^ < f -^=C qm {\\ 1 9 {r)\\ q + h 9 (r)\\ m ) 
Jo cryc — t 

(E\\ri x ' 9 (T)-rj y ' 9 (T)\\ 2 )UT < 

JO Oyjt — T Jt 

Denote by @(t) = sup X)VeG ^^^^JZ^ 3 ^'^ an d n °te that 



[j 9 {T)]a,G = SUp Xyy&G [ 



K 1 g (r)\\Vu 9 (T,x)-Vu 9 (T,y)\\ 



\x - y\\ c 
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Ku{r)\\Vg(T,x) - Vff( r, y)|| i WftW r i M r W1 



\x — y\\ a 



Then combining the above estimates for Ki, i = 1,2, 3, 4, and applying 
the Gronwall lemma we derive the estimate 

0(t) < N 9 (t)efi K tt T)dT = Nf(t), 

where N 9 (t) is a positive bounded function defined on the interval 
[0, Ti] and depending on parameters in conditions C 5.1 and C 5.2. 

The estimates of theorem 5.7 and lemma 5.8 allow to prove the 
uniform convergence on compacts of the successive approximations 
(5.10)-(5.14) for the solutions of the system (5.7) - (5.9) in C([0,Ti], 
C 1,a (K)) n C([0, Ti],L m (G) n L q (G)) forl<g<|<3<m<oo 
and arbitrary compact G in R 3 . 

To this end we differentiate the system (5.10)-(5.14) and add to 
this system the following relations 

= -Vu k (9, 4M'e dd > = J > ( 5 - 47 ) 
where / is the identity matrix acting in R? and 

Vu k+1 (t,x)=E[Vu (^\xM; k - 

Jl V +1 (r, 1% T (x)) J* Vtidw(9)dr], (5.48) 

r°° 1 

-2Vp k+1 (t,x) = -E[~f k+1 {t,x + B{T))B(T)]dT, (5.49) 

JO T 

where j k+1 = Vu k+1 Vu k . 

Now we can prove the following assertion. 

Theorem 5.9. Assume that C 5.2 holds. Then if k — > oo the 
functions u k (t),Vu)~(t,x) determined by (5.8) and (5.48) uniformly 
converge on compacts to a limiting function u(t) £ C([0, Ti], C 1,a ), < 
a < 1 /or a// t G [0,Ti], where [0, Ti] is i/te interval such that the 
solution of (5.4-5), (5.46) is bounded on [0,T\]. In addition on this 
interval the limiting function satisfies the estimates sup x \\Vu(t,x)\\ < 
a(t) , ||Vu(t)||q i i oc < (3{t) for 1 < q < | where a(t) and [3{t) solve 
(5.45), (5.46). ' 

Proof. By theorem 5.7 we know that the mapping 

«oWf,oO*0)- fvp 9 {T^L{x))dr 



$(t,x,g)=E 

acts in the space C 1,a n T g ^ oc n T m ,/oc (for a fixed i G [O.Ti]) with 
l<g<|<3<m<cx). 

Consider the successive approximations (5.10) -(5.14) and (5.47) 
- (5.49), denote by 

S k+1 (t,x) = \\u k+1 {t,x) -u k (t,x)\\, 

k+1 (t,x) = \\Vu k+1 {t,x) - Vu k (t,x)\\ 



n 

and let 



l k (t) = \\S k (t)\\oo, m k (t) = \\S k (t)\\ r , 
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p K (t) = \\n*(t)\\ 00 , Cr(t) = \\n K 
Then we obtain 

n k+1 (t,x) < Ll(E\U k tfl {x) - <o 1 WIIII<b fe ||] + 

E[\\4,o (x)\\U;} ~ <b" _1 H]) + f —^—rEiWVp^Hr, ^ >T (x))- 

Jo cr(r — rj 

v/Cr,^- 1 ^))!!!! J\lfdw(e)\\]dr+ 



1 



'T 



dr. (5.50) 



(r)dT 



o <r(t — T 

Recall that by lemmas 5.2, 5.3 we know that 
sup^H^o^-^VWII < f t [\\u k {T)-u k - 1 (T)\\ 00 ]drefo a 

x JO 

supE||<b fe - vlt'W < f HVn fc (r) - Vu k -\r)UdrJ>^ 

x JO 

+ sup f E\\Vu k -\r^l T {x)) - Vu k - l (T,4~\x))\\dreIo a ^ dT 

x Jo 

and applying the estimates from theorem 5.7 we get 

p k+ \t) < J>^ T)dT [Ll j\n V E\\u k {T^ k t:T {x))-u k -\r^ k t -\x))\\dT 



+ I" p k (r)dr + sn V f E\\Vu k -\r^ k T {x)) - V u k ~ l {r ^ k ~ l {x))\\dr] 

JO x Jo 

+ f -^=C[\\Vu k {T)Vu k -\r)\\ q + \\Vu k (r)Vu k ^(r)\\ m ] 

JO CTy/t — T 

(E\\r } k (T)- V k - 1 (r)\\ 2 00 )Ur+ 

t J t a(e)de 

/ eT r+ snpE\\Vp k+1 (r^ k T (x)) - Vp k {T^ k -\x))\\ 2 )^dT. 

JO <J\Jt — T x 

To derive the estimate for the last term we recall ( see lemma 5.1) 
that for 1 < q < | the inequality 

\\Vp k (t,x)-Vp k (t,y)\\ < ||VVWIIoolk-y|| < 

C[\h k (t)\\ q , loc + h k (t)U G ]\\x-y\\ 
holds and as a result we obtain 

E\\Vp k {r^ k T {x))-Vp k {T^ k -\x))\\ < 
C[P(r) + &(r)]EU k T (x)-^ k -\x)\\. 

In addition 

\\V P k+ \t) - V/Wlloo < c qm [\\ 1 k+ \t)- 1 k {t)\\^ loc+ 
\\l k+ \t) -VWIU.ioc] < C qm a(t)[\\Vu k+1 (t) - Vu k (t)\\ q joc+ 
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\\Vu k (t)-Vu k -\t)\\ q , loc + 
\\Vu k+1 (t) - Vu k (t)\\ m , loc + \\Vu k (t) - Vu k -\t)\\ m , loc \. 
It results from (5.50) that 

n k+1 (t,x) < C{t)[f Q E\\Vu k (T,4,r(x)) -Vu k ~\T^ k - 1 {x))\\dr+ 

f t n k (T,x)dr} + /■*_^Ci[||V« fc (r)V« fc - 1 (r)||,+ 
jo jo cry/t — t 

||V« fc (r)Vu fc - 1 (r)|| n ,] r (£;||»f' fc (r)-,f' fc - 1 (r)|| 2 )3dr 

_^ =e />(^( J B||V/ +1 (r,^ T ( : r:))-V/(r,^- 1 (x))|| 2 )^r. 
o o\Jt — t 



Note that by the Holder inequality we can prove that for any pos- 

mi r 



itive /(r) € U and -i- + \ = 1 



/ [/ /(r,x)dr] r dx< / i™i / f r (T,x)drdx = 
Jg Jo Jg jo 



i 7 ™ 1 / / f r {r,x)dxdT 
Jo Jg 



and for + £ = 1 and mi < 2 we have 



f ft | r(2-m 1 ) /■* /• 

/ [/ f(T,x)dr} r dx < t~^~ / / f( T ,x)dxdT. (5.51) 

jg jo ovt — r Jo Jg 

Then from (5.50) and (5.51) we have for r > 2 

C fc+1 W < C 2 (t)[f f [E\\u k (r,4A x ))- uk " 1 (r,^- 1 ( x Wdxdr] + 
Jo Jg 

fc, k (r)dT+ f f \\Vu k -\T^ k Ax))-Vu k ~\T^ T \x))fdxdr\ 
Jo Jo Jg 

+ f -^=C[[\\Vu k (r)Vu k - 1 (r)\\ q + ||Vu fc (r)V«*- 1 (T)IU r 

JO O-y/t — T 

( (E\\r ] x > k (T)-r l x > k - 1 (T)\\ 2 )idx}dT+ [* —t= e f' aWd0 
Jg Jo ayt — T 

[ (E\\V P k+1 (T^ k T (x)) - Vp k {T^ k - T \ X ))\\ 2 ) r *dxdT. 

Jg 

In addition for m k (t) = \\u k (t) — u k ~ l {t)\\ r j oc using the apriori 
estimates proved in lemmas 5.2 - 5.8 and theorem 5.9 we obtain 

{t)<C{t)[{f I E\\u k (r^ k T (x))-u k - 1 (r^ k -\x)Wdxdr)T 

J G 

+(1^-5 f [ E\\Vu k +\T^ k +\x))Vu k {T^ k T {x))- 
o Jo Jg 

Vu k (r^ k T (x))Vu k -\T,ij k - 1 (x)WdxdT) 1 -]< 
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m k+1 



Ci(t)[( /' m k r (T)dT); + ( /' / a(r)S||^ T (x)-^ 1 ( a; )ir da;dr )' + 
jo jo Jg 

-t^{[\p k+ \T)+p k {T)]^{T)dT)-r]. 

C Jo 

Since u k and Vu fc are uniformly bounded on [0, Ti] and 
||Vii 1 (t, •) — Vuo(-)\\ r ,ioc < const < oo, 

both for r = m and r = q we obtain that there exists a bounded on 
[0, Ti] positive function C2(t) such that the function K n (t) = p n {t) + 
Cm(t) + m r satisfies the estimate 

«»«) < |C2W1 " 



n! 

and hence lim^oo «; n (t) = 0, since C^i) is bounded on [0,Ti]. Fi- 
nally we obtain that for each i G [0, Ti) the family it n (t, •) uniformly 
converges to a limiting function •) G C 1,Q nL m j oc . In addition, we 
can check that the limiting function Vu(i, x) is Lipschitz continuous 
in x. In fact, by lemma 2.8 and theorem 2.9 for each t G [0,Ti], we 
have for any x,y £ G 

\\Vu n (t,x) -Vu n (t,y)\\ < N(t)\\x - y\\, 

where iV(t) and Ti were defined above in lemmas 5.8 and theorem 5.7 
respectively and the estimate is uniform in n. This allows to state 
that the limiting function is Lipschitz continuous as well. 

To prove the uniqueness of the solution of (2.8)-(2.10) constructed 
above we assume first that there exist two solutions ui(t, x), U2(t, x) to 
(5.7)-(5.9) possessing the same initial data «i(0, x) = «2(0, x) = uq(x). 

Computations similar to those used to prove the convergence of 
the family (u n (t),Vu n (t)) allow to check that 

(Vui(t)-V«2(t)]a,G = and || V«i(t) - Vu 2 (t)\\ m ,i oc = 0. 

Finally, we know that the Cauchy problemfor a stochastic equation 
with Lipschitz coefficients has a unique solution. This implies the 
uniqueness of the solution to (5.7)-(5.9). 

Summarizing the above results we see that the following statement 
is valid. 

Theorem 5.10. Assume that C 5.2 holds. Then there exists 
a unique solution ■ip tjX (s),u(t,x),p(t,x) of the system (5.7)-(5.9), for 
all t from the interval the [0, Ti], with T\ given by theorem 5.7 and 
x G G for any compact G C R 3 . In addition ipt,x(s) is a Markov 
process in R 3 and u G C([0,Ti], C 1 ' a (G)) n C([0, Ti], L qM n L m>loc ) 
for l<q<|<3<m<oo. 

To fulfill our program we have to check that the conditions of the- 
orem 2.8 are sufficient to verify that the functions u(t, x),p(t, x) given 
by (5.8), (5.9) define a weak solution of the Navier -Stokes system. 

To this end we have to apply the results of the Kunita theory of 
stochastic flows. 
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Namely we check that given a distribution valued process A(i) of 
the form 



A(t) = no - / Vp u (r) o $ T dr (5.52) 

JO 



the function 



A(t) o = o iPtfl - J* Vp u {r) o ^ T dr 

gives rise to a solution of (5.1). 

To this end we apply the generalized Ito formula [19], [20] to derive 

A(t) o = + f Q yAKO) o ^ o ]d0+ (5.53) 

fv[u{9) o^ Q ]adw(d) - f t V[u(9)o^ ] u (e)de- f Vp u {e)d9. 
Jo Jo Jo 

2 

Note that for Lu = — (u, V)u + ^-Au we have 



JO 



E 

Hence 



(L(u(t) o V" s (x))dr, /i(a;))dx = 
jr l <u(r)o^ )0 ,L*/»)dr] = j\{E[u{ror T ,o)lh)dT. 



u(t) = E[X(t) o = «o + f LE[u(t) o ^« ]dr - /' Vp"(r)dr. 

JO JO 

Differentiating each term with respect to t we can check that the 
function 

u(t)=E[\(t)o^ ] (5.54) 

solves the Cauchy problem (5.1), (5. 2). 

To summarize the obtained results we can state the following as- 
sertion. 

Theorem 5.11. Assume that C 5.2 holds. Then the functions 
u(t,x), p(t,x) given by (5. 8), (5. 9) are defined on the interval [0, T\] 
with T\ determined by theorem 5.8 and satisfy (5.1)- (5. 2) in a weak 
sense on this interval. 

Remark 5.12. We have proved that under condition C 5.2 the 
system (5.7)-(5.9) gives rise to a weak solution of (5.1)-(5.2) . Note 
that if the initial data are smoother, say uq £ C 2+a , a £ [0, 1] similar 
considerations can be applied to verify that the pair u(t,x),p(t,x) 
given by (5.8)-(5.9) stands for a classical C 2 -smooth solution of (5.1)- 
(5.2). In fact in this case applying the generalized Ito formula for the 
verification assertion we may treat the action of the operator L in the 
classical sense rather then in the weak sense. 
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6 Lagrangian and stochastic approach 
to the Euler and the N-S system 

The probabilistic approach developed in the previous section is in 
a sense an analogue of the Lagrangian approach to the Euler and 
the Navier-Stokes systems. A rather close model was constructed in 
papers by Constantin and Iyer [11], [12]. To make it easier to compare 
we rewrite the results from these papers in terms similar to those used 
in the previous section. We consider first the Euler system 

^ + (u, V)u = -Vp, u(0, x) = u (x), x € R 3 (6.1) 

divu = 0. (6.2) 

and recall that the corresponding Lagrangian path starting at y is 
governed by the Newton equation 

%^ = ^). (6.3) 
The incompressibility condition for the map <j) yields 

det(V<My)) = 1. (6.4) 
The force F in (3.3) has the form 

Ffo, y ) = -v P (t,4>oAy)) = -[(VMvVT^WiMy))]- ( 6 - 5 ) 

One can deduce from (6.3) that 

d d% tt (y)d%, t (y) dgjtMM) (RR , 

dt [ dt d Vl 1 d Vi ' 1 j 

where 

,(t,y)=P(t,y) (6.7) 

We recall that in (6.6) and below summation over the repeated indices 
is assumed. Integrating (6.6) in time we get 

d~^ t {y)d^ t {y) _ dn&Mv)) ( , R) 

~di dy-- Uoiy) % — ' (6 - 8) 

where ^ 

n(t,y) = I q(r,y)dT (6.9) 
Jo 



My) = — t — \t=o (6-io) 



and 

= — - 

dt 

is the initial velocity. 

Consider the inverse diffeomorphism ipt,o = [<^o,t] _1 come back to 
(6.7), multiply it by [V^,o] and put y = ipt,o(x)- As a result we obtain 
by the chain rule the relation 

u% x) = (v&$t fi (x))V Xi $ fi (x) - f V Xi q(r, M x )) dT - (6-H) 

Jo 
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The equation (6.11) shows that the general Euler velocity may be 
written in the form that generalizes the Clebsch variable representa- 
tion 

u = [VVvjfC - Vn, 

where C = uo(tptfi(x)) is an active vector and n is defined by the 
incompressibility condition divu = 0. 

Note that a vector A is called active if 

-A = - + (u,V)A = . 

It is easy to check by the chain rule that 

J t M*) = + («-V)^W = 0, (6.12) 

that is vpt,o(x) is an active vector. 

Hence the Euler equations are equivalent to the system consisting 
of (3.9) and the following relation 

An(t, x) = ±{4(A :0 (x))^^}, (6.13) 

where n is given by (6.9). 

Now one can assume the periodic boundary conditions or the 
zero boundary conditions at infinity. Note that in the periodic case 
n(t, x),u(t, x) and 

8{t,x) =x-4> tfi {x) (6.14) 

are periodic functions in each spatial direction. Finally due to divu = 
one can rewrite the equation of state (6.11) in the form 

u{t) = n{n J (^,o)V^ } = n{[V^ l0 ]*«o(^,o)}, (6-15) 

where II = I — VA _1 V is the Leray-Hodge projector (with corre- 
sponding boundary conditions) on divergence free vector fields. The 
Euler pressure is determined up to additive constants by 

p(t,x) = dU ^ x) + (u(t,x),V)n(t,x) + ^\\u{t,x)f. 

Note that (6.11), (6.12) made a closed system and may be used to 
determine u{t). 

Let us compare (6.11), (6.12) with the alternative representation 
for the state u{t) developed in the previous section. 

To this end we choose 4>o,t '■ U — > 4>o,t{y) to be a volume preserving 
diffeomorphism that satisfies the equation 

#o,r(y) = u(t - r, 4> 0:T (y))dT, 4>o,o(y) = y, (6.16) 

with div u(t) = 0. 

Consider the system 

dipt,d(x) = -u(9,ip t o{x))d9, ip t ,t( x ) = x > ( 6 - 17 ) 
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u(t,x) = u (i>t,o(x)) - [ Vp(T,ipt,A x )) dT i ( 6 - 18 ) 

Jo 

roo 

-2p(t,x)=E[ ~f{t,x + B{T))dT], (6.19) 
Jo 

where 7 is given by (1.3). 

If the fields u(t,x),p(t,x) are regular enough then we may con- 
struct the representation of the solution to the Euler system in the 
form (6.18), (6.19). 

To check this we consider a volume preserving diffeomorphism </>o,t : 
V </>o,t(y) that satisfies (6.16). 

Next we consider the system 

dipt,e(x) = -u(6,ip tj o(x))d9, ip t ,t{ x ) = x > ( 6 - 20 ) 

u(t,x) = u {ip t ,o( x )) - [ ^Pi T AtA x )) dT ^ ( 6 - 21 ) 

Jo 

roo 

-2p(t,x) = E[ j(t,x + B(T))dr], (6.22) 
Jo 

where 7 is given by (1.4). 

If u(t, x) is regular enough then we may construct the representa- 
tion of the solution to the Euler system (6.1), (6.2) in the form (6.21), 
(6.22). 

To this end we consider a vector field \{t) satisfying the equation 
^1 = -Vp(t) o O;t A(0) = «o, 

where <pQ t satisfies the ODE (6.16), and let the process ipt,o be its 
inverse. Applying the Kunita approach [19] to the process ip t ,o we can 
verify that Y't.o along with (6.20) satisfies the equation 

Yv(z) = x + j* V^ )t (^ t)0 )- 1 u(^, x)d6, (6.23) 

that allows to prove that u(t) given by (6.21) satisfies (6.1). 

Comparing (6.15) and (6.21) we note that they give different ex- 
pressions for the velocity field. Actually (6.21) includes the Euler 
pressure p(t,x) instead of q(t,x) used in (6.15). Besides the proba- 
bilistic representation of the solution p to the Poisson equation 

-Ap = ViU k V k Ui 

is used instead of the Leray projection. 

Coming back to the Navier-Stokes system ((5.1), (5. 2) we recall 
here the approach due to Constantin and Iyer [12]. The stochastic 
counterpart of the Navier-Stokes equations in the version of Iyer [25] 
looks like the following. 

Consider the closed stochastic system 

#0,0 = u{9, <j>Qfi)dt + adw(9), (/) ,o{y) = y, (6.24) 
^,0 = [<M~\ (6.25) 
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u(t) = M[(V^ t ,o)(tio o ^ j0 )]. (6.26) 



The existence and uniqueness of the solution to this system is proved 
in [12] by the successive approximation technique. As a result the 
authors constructed a strong local in time solution of the Cauchy 
problem for the Navier-Stokes system for regular enough initial data. 

The main result due to Constantin and Iyer reads as follows 

Theorem 6.1. Let k > 1 and uq G C k+1 ' a be divergence free. 
Then there exists a time interval [0, T] with T = T(k,a,L, ||uo||fc+i,a) 
but independent of viscosity a and a pair <^o,j(^), u(t, x) such that u £ 
C([0,T],C fc+1 ' Q ) and (u,<f>) satisfy (6.24)-(6.26). Further there exists 
U = U(k, a, L, \\uo\\k+i : a) such that \\u(t)\\k+i t a < U for t £ [0, T] and 
u satisfies the N-S system. 

As we have mentioned above an approach close to the one of [12] 
was developed in our previous paper [7]. Both these approaches allow 
to construct a classical (local in time) solution to the Cauchy problem 
for the Navier-Stokes system and prove the uniqueness of the solution. 

On the other hand the approach developed in section 5 allows to 
construct a weak (local in time) solution to (5.1), (5.2) and prove the 
uniqueness of this solution in the corresponding functional classes. 

The stochastic counterpart of the Navier-Stokes system considered 
in section 2 has the form 

dip t ,e(x) = -u(9, ip tfi (x)d9 + adw{6), ipt,t( x ) = x > (6.27) 

u(t, x) = E[uo(ilH fi (x)) - f Vp(r, ik tT (x))dT], (6.28) 

J o 

roo 

2p(t, x) = - E[y(t, x + B(r))}dT. (6.29) 
Jo 

Note that we can use the relation 

-2Vp(t, x) = E[ / -j(t,x + B(T))B(T)dT] (6.30) 

JO T 

to eliminate the pressure from the above system (6.27) - (6.29). 

We can see that the difference between (6.27) - (6.29) and (6.24) 
- (6.26) has the same nature as the difference between (6.11), (6.12) 
and (6.17) - (6.19). 

Finally we note that the approach developed in section 5 allows us 
to construct both strong (classical) and weak (distributional) solutions 
of the Cauchy problem for the N-S system. 

Acknowledgement. The authors gratefully acknowledge the fi- 
nancial support of DFG Grant 436 RUS 113/823. 

References 

[1] McKean H. A class of Markov processes associated with nolinear 
parabolic equations. Proc. Nat. Acad. Sci. USA 59 6 (1966) 1907- 
1911. 



46 



[2] Freidlin M. Quasilinear parabolic equations and measures in func- 
tional spaces. Funct. Anal, and Appl. 1, N 3 (1967) 237-240. 

[3] Freidlin M. Functional integration and partial differential equa- 
tions. Princeton Univ. Press 1985. 

[4] Belopolskaya Ya., Dalecky Yu. Investigation of the Cauchy problem 
for systems of quasilinear equations via Markov processes . Izv 
VUZ Matematika. N 12 (1978) 6-17. 

[5] Belopolskaya Ya.I. Dalecky Yu.L. Stochastic equations and differ- 
ential geometry. Kluwer (1990). 

[6] Kunita H. Stochastic flows and stochastic differential equations, 
Cambridge Univ. Press, Cambridge. (1990). 

[7] Albeverio S., Belopolskaya Ya. Probabilistic approach to hydrody- 
namic equations. In the book "Probabilistic Methods in Hydrody- 
namics." World Scientific (2003) 1-21. 

[8] Elworthy K. D., Xue-Mei Li, Differentiation of heat semigroups 
and applications, J. Funct. Anal. 125, No.l, (1994) 252-286. 

[9] Albeverio S., Belopolskaya Ya., Probabilistic interpretation of the 
VV- method for PDE systems. In the book "Analytical approaches 
to multidimensional balance laws", Nova Science Publishers, Inc. 
(2005) 1-43. 

[10] Albeverio S., Belopolskaya Ya., Probabilistic approach to systems 
of nonlinear PDEs and vanishing viscosity method Markov pro- 
cesses and related topics 12 1 (2006) 59-94. 

[11] Constantin P., An Eulerian-Lagrangian approach to the Navier- 
Stokes equations. Commun. Math. Phys. 216 (2001) 663-686. 

[12] Constantin P., Iyer G. A stochastic Lagrangian representa- 
tion of the 3-dimensional incompressible Navier-Stokes equations. 
arXiv:math. PR/051067 4 31 Aug 2006. 

[13] Leray J. Essay sur le movement d'un fluide visqeux emplissant 
I'espace, Acta Math. 63 (1934) 193-248. 

[14] Kato T., Strong L p solutions of the Navier-Stokes equations in 
R m with application to weak solutions. Math. Zeit. 187 (1984) 
471-480. 

[15] Lemarie-Rieusset P.G. Recent Developments in the NavierStokes 
Problem. Chapman&Hall / CRC (2002). 

[16] Le Jan Y., Sznitman A.: Stochastic cascades and 3-dimensional 
Navier Stokes equations, Prob. Theory Relat. Fields, 109 (1997) 
343-366. 

[17] M. Ossiander A probabilistic representation of solution of the in- 
compressible Navier-Stokes equations in R 3 . Prob. Theory Relat. 
Fields, 133 2 (2005) 267-298. 

[18] B. Busnello, F. Flandoli, M. Romito, A probabilistic representa- 
tion for the vorticity of a 3D viscous fluid and for general systems 
of parabolic equations, Proc. Edinburgh Math. Soc, 48 2 (2005) 
295-336. 



47 



[19] Kunita H. Stochastic flows acting on Schwartz distributions. J. 
Theor. Pobab.7 2 (1994) 247-278. 

[20] Kunita H. Generalized solutions of stochastic partial differential 
equations. J. Theor. Pobab.7, 2 (1994) 279-308. 

[21] Belopolskaya Ya., Woyczynski W. Generalized solutions of non- 
linear parabolic equations and diffusion processes. Acta Applican- 
dae Mathematicae 96 1-3 (2007) 55-69. 

[22] D.Gilbarg, N.S.Trudinger. Elliptical partial differential equa- 
tions second eddition Springer- Verlag Berlin Heidelberg New York 
Tokyo (1983). 

[23] Ikeda N., Watanabe S. Stochastic differential equations and dif- 
fusion processes, North-Holland Publ. Comp. (1989) . 

[24] Constantin P., An Eulerian-Lagrangian approach for incompress- 
ible fluids: local theory. JAMS 14 (2001) 263-278. 

[25] Iyer G. A stochastic Lagrangian formulation of the incompressible 
Navier-Stokes and related transport equations. PhD dissertation 
Dept. Math. Univ. Chicago (2006). 



48 



